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ABSTRACT 

Context. Since the external regions of the envelopes of rapidly rotating early-type stars are unstable to convection, a coupling may 
exist between the convection and the internal rotation. 

Aims. We explore what can be learned from spectroscopic and interferometric observations about the properties of the rotation law in 
the external layers of these objects. 

Metliods. Using simple relations between the entropy and specific rotational quantities, some of which are found to be efficient at 
accounting for the solar differential rotation in the convective region, we derived analytical solutions that represent possible differential 
rotations in the envelope of early-type stars. A surface latitudinal differential rotation may not only be an external imprint of the inner 
rotation, but induces changes in the stellar geometry, the gravitational darkening, the aspect of spectral line profiles, and the emitted 
spectral energy distribution. 

Results. By studying the equation of the surface of stars with non-conservative rotation laws, we conclude that objects undergo 
geometrical deformations that are a function of the latitudinal differential rotation able to be scrutinized both spectroscopically and by 
interferometry. The combination of Fourier analysis of spectral lines with model atmospheres provides independent estimates of the 
surface latitudinal differential rotation and the inclination angle. Models of stars at different evolutionary stages rotating with internal 
conservative rotation laws were calculated to show that the Roche approximation can be safely used to account for the gravitational 
potential. The surface temperature gradient in rapid rotators induce an acceleration to the surface angular velocity. Although a non- 
zero differential rotation parameter may indicate that the rotation is neither rigid nor shellular underneath the stellar surface, still 
further information, perhaps non-radial pulsations, is needed to determine its characteristics as a function of depth. 

Key words. Stars: early-type; Stars: rotation; Stars: spectroscopy; Stars: interferometry 

1. Introduction Nevertheless, many attempts have been made to obtain in- 

formation on the internal rotation from detailed studies of: a) the 
1.1. Review of observational approaches position of stars in the HR diagram; b) the evolution of the Vsin i 

One of the most enduring unknowns in stellar physics has been parameter during the main sequence (MS) phase; c) the shape of 

" " the inner distribution of the angular momentum in a star. In absorption lines, whose characteristics can depend upon the ro- 

the past few decades, significant progress has been made in de- ^^^10"^! in layers close to the stellar surface; d) the global 

scribing theoretically the evolution of rotating stars. This has re- stellar geometry described with mterferometnc data, 
quired an understanding of numerous hydrodynamic and mag- briefly leview these effoits. 

netic instabilities triggered by the rotation, as well as the mixing a) The most numerous eff^orts among those just mentioned 

processes of chemical elements unleashed by these instabilities are the statistical analysis of photometric data on die rotational 

fesoul 1978 .2000; Zahii. 1983 , , 1992; Maeder & Eenens 2 004^ spread of the MS, which can be described bv (iRoxb urgh et al 



lMaedeiil2009h . However, apart from flie Sun, reliable observa- 1956; Maedd[T%l ICollins & SonnebornlligTTt ICollins et^ 



tional information about the internal rotation of stars remains ( 199 id 
scarce or non-existent. 

AMv = k(n)V", (1) 

Send offprint requests to: J. Zorec: e-mail: zorecOiap . fr 

* Member of the Carrera del Investigador Cienti'fico, CONICET, where AMy is the deviation in absolute magnitude from the 

Argentina zero-rotation MS, V is the true equatorial rotational velocity. 
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and k(n) is a constant whose value depends on the power n. 
When n - 2, k{n = 2) is on the order of kglO^^ mag/(km s"')^, 
so that for k„ < 1 the deviations may indicate that the in- 
ternal rotation is uniform, while for k„ > 1 the internal ro- 
tation can be differential dCotton & SmithI Il983h . This type 
of analysis found that stars do not seem to rotate uniformly. 
However, owing to the measurement uncertainties and dif- 
ficulties in defining the MS of zero rotation, the available 
data could not provide a ny firm evidence of a particu l ar law 



aata couia not proviae a ny nrm evidence or a particu l ar law 
of non-uniform rotation ('Strittmatter & Sargent 119661: Golayl 
1968 ; Maeder 1968; Maeder & Peytremann||197a |SmM]|[192lF 



Smith & Worlevlll974t iMoss & Smith 1981). Furthermore, us- 
ing detailed model atmo spheres for differentially rotating stars, 
ICoUins & SmithI ( Il985h concluded that photometry alone can 
place albeit rather weak constraints on the degree of differential 
rotation within the stars. 

b) Depending upon the internal angular momentum redistri- 
bution and evolutionary rearrangements of the inertial momen- 
tum, the surface equatorial rotational velocity of stars changes 
accordingly. Thus, the study of the variation in the true rotational 
velocity, V, as a function of time was studied by several authors 
using the ratio 



(Vsin />Lc 



<V>LC 



(Vsin />ZAMS < V)zAMS 



(2) 



where (Vsin /)lc is the average of the Vsin / parameters of stars 
with in principle the same mass and luminosity class (LC), 
(ysini)zAMS is the average of Vsin / for stars with the same 
mass, but located near the zero-age-main sequence (ZAMS). 
These ratios were compared with similar ones predicted theo- 
retically for stars evolving as rotators in two different and ex- 
treme ways. On the one hand, the stars were assumed to evolve 
all their way as uniform rotators, which implies that the angu- 
lar momentum is entirely redistributed at each evolutionary step. 
On the other hand, it was assumed that each stellar layer con- 
served its initial specific angular momentum, i.e., the stars did 
not undergo any redistribution of its internal angular momen- 
tum. Since in many cases the observed ratios 7?lc were found 
to be situated in-between the two extreme theoretical predic- 
tions, it was suggested that stars should be diff'erential rota- 
tors. However, these studies could not provide any information 
about the characteristics of the internal rotational law ( Sandaed 



aoout tne cnaracteristics or tne internal rotational law (laanaagel 
Il955t iDanziger & FabedfT97l Iz orec et al ] ll987t IZoredl2004 . 
Somewhat related to this category of inquiries is the study of 
the evolut ion of the to t al ang ular momentum of B and Be stars 
carried bv lZorec et al.l (Il990l) . who concluded that these objects 
should undergo some internal angular momentum redistribution 
to explain the observed evolution of the Vsin ; parameters. 

Ci) The study of the absorption line profiles of MS B-type 
stars found some evidence of possible surface diff'erential rota- 
tion. The angular velocity in the surface of stars was assumed to 
depend on the colatitude angle, 6, as 



Q(6») = QJl - 5 X (1 - R(0) sin^ ff)] 



(3) 



where Q.o is the equatorial angular velocity, R{0) is the equation 
of the stellar surface, S is the parameter that testifies to the differ- 
ential rot ation. Using stellar mod els that are more or less gravity 
darkened. lStoecklevI (1 1 968ah and lStoecklev & Bu scombd (ll987l) 
found that in most cases S < 0, which suggested that the an- 
gular velocity tends to increase from the equator to the pole. 
Nevertheless, a dependence of the surface angular velocity on 



the latitude Q(0) could be due either to an actual differential 
rotation present under the stellar surface, or simply to zonal 
atmospheric currents, which cou ld appear in rapidly rota ting 
early-type stars, as speculated bv iCranmer & CollinsI (Il993h . In 
Sect. 12.41 we recall that an acceleration of the angular velocity 
towards the equator, i.e. 5 > 0, can be promoted by a temperature 
gradient induced by the gravity darkening effect. 



C2) The possibility of detecting surface differential rota- 
tion by means of the Fourier a n alysis of spectral line profile s 
was discussed by Huang (1961), Gray ( 1977), Bruning (19811), 
Garcia- Alegre et al. (1982), and Reiners & Schmitt (200J. 
Evidence of surface differential ro tation in late-type s t ars wit h 
Vsi n i < 50 km s were given by IReiners & Schmitt" ('2003bl), 
Rei ners & Schmitt! (|2004 . and joiners & Royer (2004), but no 
differential rotation for late-type stars with Vsin / > 50 km s"' 
and A-type stars with Vsi n / > 150 km s~' were reported by 
Rei ners & SchmittI (l2003al) and iGravl (1 19771) . respectively. It is 
possible that modest differential rotation is difficult to detect 
with the Fourier transform technique in slowly rotating A-type 
stars, because the rotational broadening is not large compared 
with the broadening caused by other mechanisms such as ther- 
mal turbulence and pressure effects riravlll977l) . However, in 
those cases where there is some evidence of differential rota- 
tion, the parameter S cannot be differentiated from the unknown 
inclination angle factor sin /. Nevertheless, its sign indicates ac- 
celeration of the angular velocity towards the equator. To our 
knowledge, the Fourier technique for differential rotation has not 
yet been applied to early-type stars. 



C3) lAndd (Il980l) proposed a method to probe the inner an- 
gular velocity of stars based on the use of the rotational split- 
ting of non-radial oscillations. However, owing to uncertainties 
in the identification of pulsation modes and rough determina- 
tions of stellar fundamental parameters, particularly their evolu- 
tionary stage, this method has not yet been able to be applied 
with reliable success. Nevertheless, from the analys is of pulsa- 
tion modes derived from photometric variations ( Dupr et et al.l 
2004) and COROT data (Degroote et al. 2009), constraints on 
the internal rotation of /3 Cep stars have been inferred. 



d) In the past few years, interferometric methods have helped 
provide remar kable insights into not only the rotational distor- 
tion o f stars (iDomiciano de Souza et al.l l2003t Ivan Belle et al.1 
20041) . but also the induced gravity darkening effe ct by means 
of imaging technique s (McAlisteret al. 2005; Aufd enberg et al.l 
2006; van Bel le et al.l l2006; Zha o et alJl200 9). New instruments 
with higher spectral resolutions of up to the 10000 attained by 
VLTI/ AMBER in the J and K b ands and an angul ar resolution 
of about 1 mas in the K band ( Petrov et alJl2007l) . or spectral 
resolution reaching 30000 and angular resolutions as high as 0.3 
mas in the visible using the VEGA/CHARA interferometric ar- 
ray (Mourard et al. 2009). will not only probably enable us to 
determine with greater detail than in previous studies the global 
geometry of stars deformed and gravity darkened by the rotation, 
but also carry out differential interferometry. 

A method based on differential interferometry that re- 
quires high spectr al and spatial resolution was presented by 
Domiciano de So uza et al.l y004b c a) to distinguish observa- 
tionally the parameter controlled by the degree of the surface 
differential rotation from the inclination angle factor sin 
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1.2. Aims of the present attempt 

Most theoretical predictions about the evolution of rotating 
early-type stars and the mixing of chemical species triggered 
by the instabilities set up by the rotation, come from calcula- 
tions performed in the framework of two significant assump- 
tions: a) the global rotational energy stored by the stars in the 
ZAMS is lower than the limit allowed by the rigid rotation in the 
critical regime; b) the internal angular velocity undergoes an in- 
stantaneo us "shellular" re distribution at each evolutionary step. 
However, iciemen ll (fl979l) using a c yhndrical (conservati ve) dis- 
tribution of the angular velocity, and'Mae der et alJ (l2008h basing 
their calculation on a "shellular" distribution, showed in a more 
detailed way that in rapidly rotating early-type stars the envelope 
layers beneath the surface, may have wider convective zones 
in radius than in non-rotating stars. In the Sun, only the layers 
unstable to convection rotate differentially with a non-shellular 
pattern. This motivates the inquiry of whether in massive and 
intermediate mass stars some coupling may also exist between 
convection and rotation beneath their surface. In this case, the 
characteristics of the rotational law in the external stellar layers 
should differ from those currently assumed in the above evoked 
stellar models. 

As demonstrated by many authors, the global geometry of 
a star depends not only on the total amount of angular mo- 
ment um stored by the s t ar, but also on its internal distribu- 
tion ('Bodenheimer"l97l'; 'Zored ll986t ISmith & CoU"inslll992t 
[Urvu & Eriguchi 1994, 1995). This geometry mostly relies on 
the stellar surface rotation, which acts as an imprint of its prop- 
erties in the layers beneath the surface. In this case, we should 
not exclude the resulting mixing of chemical elements in the stel- 
lar atmosphere being more or less dependent on the characteris- 
tics of the external rotational law, upon which the description 
of the stellar structure, based on the abundance determination 
of chemical elements, should also rely. Therefore, to provide 
new information and/or constraints to test the global assump- 
tions currently made to calculate models of stellar structure with 
rotation and thus help deepen our understanding of the proper- 
ties of early-type fast rotators, we might ask: 1) what can be 
deduced, using first principles, about the properties the rotation 
laws can have beneath the surface as a consequence of the cou- 
pling between rotation and convection; 2) what are the parame- 
ters needed to characterize these stars that may be accessible to 
observations; 3) whether the combined interpretations of spec- 
troscopic and interferometric data of rapidly rotating early-type 
stars enable us to determine these parameters. 

In this attempt, the most interesting information might be the 
indication of some differential rotation in the stellar surfaces and 
the sign of its latitudinal gradient. Both pieces of data can be ob- 
tained, as much as possible, in a consistent way by taking into 
account the stellar geometrical deformation produced by this ro- 
tation and the concomitant gravitational darkening effect that re- 
sponds to possible non-conservative rotation laws. 

The present paper is organized as follows. In Sect. |2] we use 
first principles to infer possible rotation laws in the convective 
layers beneath the stellar surface of early-type rapid rotators. 
Sect. [3] presents the equation of the stellar surface of stars with 
non-conservative rotation laws. A discussion of the gravity dark- 
ening effect for non-conservative rotation laws is presented in 
Sect. |4] The discussion about the validity of the Roche approx- 
imation in representing the gravitational potential is presented 
in Sect. |5] This discussion is based on 2D models of rotating 
stars where the evolutionary stages are taken into account in a 
simplified way. We briefly comment on the determination of the 



rotational profile in the stellar envelope in Sect.|2l In Sect.|6l we 
summarize the attainable information on rapidly rotating early- 
type stars with external differential rotation from the combined 
analysis of spectroscopic and interferometric data. Our conclu- 
sions are presented in Sect. |8] 

2. Rotational law in the stellar envelope 

2.1. The angular velocity distribution beneath the stellar 
surface 

The effects of rotation are generally introduced in the struc- 
ture equations of rotating stars by replacing the spherical strat- 
ification of non-rotating star-models by a rotationally distorted 
stratification, which keeps the whole calculation problem to one 
dimension (Kippenhahn & Thomas 1970; Endal & Sofia 197a 

Pinsonneault et al. II990I: Flieaner & Langej 

,1995: ,Chabover et aDflggllMevnet & ^faeder,1997) . This pro- 
cedure is justified if the internal differential rotation has a shellu- 
lar distribution law because it comes from theoretical inferences 
made by Zahn (1983, 1992), which rely on the assumption that 
the horizontal turbulence is much stronger than the vertical one. 
The limiting case of a shellul ar rotational profi le is rigid rota- 
tion. From the calculations bv lMaeder & Mevn et (2005), it fol- 
lows that magnetic fields created by the Pitts & Tayler instability 
(iPitts & Ta vler 1985; Spruit 1999, 2002) can lock the stellar lay- 
ers to eac h other and force the star to evolve as a rigid rotator. 
However, IZahn et aTl ( l2007h concluded that the dynamo action 
can be less efficient as previously expected (Spruit 1999, 2003) 
and that the magnetic fields created contribute little or less to the 
transport of the angular momentum. In addition, if the magnetic 
field in early-type rapid rotators is created, it could perhaps has 
some effect on their rotational profile in the convective parts of 
the envelope as in the Sun (Balbus 2009; Balbus et al. 2009). 

The differential rotation in the surface of the Sun is a direct 
consequence of the differential rotation in the convective layers 
beneath the surface layers. Therefore, it is important to stress that 
this differential rotation in depth is not shellular, in spite of the 
strong turbulent viscosity that according to the above-mentioned 
theoretical assumptions, might otherwise cause shellular-like ro- 
tation^ 

iMaeder et al.l (l2008l) demonstrated that rotation does not in- 
hibit convection as could be thought from Solberg-H0iland's sta- 
bility criterion, but it changes the thermal gradient so as to en- 
hance convection. Hence, in rapidly rotating massive stars the 
two external convection zones, associated with increased opaci- 
ties due to He- and Fe-ionization, respectively, are both consid- 
erably enlarged in depth so that the entire convective zone covers 
a non-negligible external region, which ranges from some 1/8 of 
the stellar radius in the pole to nearly 1/4 in the equator. 

However, the angular momentum distribution in the convec- 
tive regions remains a puzzling question. To account for it in 
the stellar model calculations, two extreme approximations have 
been used: a) rigid rotation, angular velocity Q. = c onstant, which 
is supposedly promoted by the turbulent viscosity (iMaeder et alJ 
2008); b) constant specific angular momentum j (J - Qw^; ru 
is the distance to the axis of rotation), possibly due to the re - 
distribution produced by the convective plumes (lTavleijll973h . 
Nevertheless, the solar convective regions, characterized by sig- 
nificant turbul ence, are the on ly ones with significant differen- 
tial rotation ( Scho u et a l]fT998h . even though the solar rotational 
profile in the convective region, revealed by the helioseismolog- 
ical data, does not fall between these two extreme possibilities 
(lDeupreell200lh . Global insights into the rotational law in the 
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12.08 12.1 12.12 12.14 12.16 12.18 12.2 12.22 12.24 

log R (cm) 

Fig. 1. Adiabatic gradient d In Pp^^ ld\n r in the two rotationally 
enlarged convective regions in the envelope of a 20Mq star. The 
gradients are shown as a function of the logarithm of the stellar 
radius in the polar and equatorial directions. 



convective regions could be obtained by exploring the solutions 
to the thermal wind balance relation under imposed conditions 
between the entropy and rotation. This procedure succ essfully 
explains the Solar rotatio n law in the convective regions (iBalbusI 
l2009HBalbus et al.112009 '). However, we pay attention here to our 
dealing with rapidly rotating stars. Detailed physical justifica- 
tions of the assumptions made in the present paper are postponed 
to future contributions. 



2.2. Rotation inferred from tiie baroclinic balance relation 

The curl of the time-independent momentum equation of an 
inviscid, axisymmetric rotating star with negligible magnetic 
fields, yields the baroclinic balance relation. Using the cylindri- 
cal co ordinates (m, (p, z), this balance condition reads dTassoull 
[T978h 



= ^(VPxVp).e^, 



(4) 



where j is the specific angular momentum (j - Qtu^), P and 
p are pressure and density, respectively, and e^ is the azimuthal 
unit vector. Taking into account the equation of hydrostatic equi- 
librium 



-VP = geff, 
P 



(5) 



where geff is the effective gravity, and using for the specific en- 
tropy S the expression 



r-1 



InPp + constant. 



(6) 



where k is the Boltzmann constant and y is the ratio of specific 
heats at constant pressure and constant volume per unit mass, the 
wind equation in Eq. (|4]i can be rewritten as 



1 df 1 
— ^ = 7^(V5 xgeffO.e^, 
OZ Cp 

where Cp is the constant pressure specific heat. 



(7) 



We can attempt a discussion of Eq. (|4]i by seeking solutions 
for the stellar internal rotation under at least three different con- 
ditions where for the moment the effects carried by the merid- 
ional circulation are neglected: a) marginal stability imposed 
by the Solberg-H0iland criterion; b) state enforced by parallel 
surfaces of specific entropy and specific angular momentum; c) 
frame where the specific entropy parallels the local specific ki- 
netic rotational energy. 

a) The Solberg-H0iland stability criterion states that "a baro- 
clinic star in permanent rotation is dynamically stable against ax- 
isymmetric perturbations if two conditions are satisfied: (i) the 
specific entropy S increases outwards, and (ii) on each surface S 
- constant, the specific angular momentum increases from the 
pole to th e equator". Th e second condition is written mathemat- 
ically as (iTassoulll 19781) 



-^^[V/xV5] >0, 



(8) 



where g. is the z-component of the effective gravity. Making the 
ansatz for strict marginal equilibrium, we have 



(9) 



Equation (|9]l suggests then that surfaces = constant and 
S - constant should be parallel, S = S{j^), i.e. a displaced 
fluid element in baroclinic turbulence retains both entropy and 
angular momentum. We do not consider here the conditions that 
might render possible the balance implied by Eq. (|9]l. 

b) The assumption that the surfaces of constant specific en- 
tropy and constant angular velocity coincide, i.e. S - S(Q.^), 
brings another alternative solution to the baroclinic equilibrium 
equation (|4]i. This coupling can be enforced by magnetic fields 
(lBalbusll2009.) . altho ugh hydrodynamic constraints in the Sun 
can justify it entirely (iBalbus et al.ll2009h . 

c) Only as an extrapolation to the conditions S - S (f) and 
S - S {Q?) may we also consider S - S (nr^il^), since two energy- 
related quantities are parallel: specific entropy S and specific ro- 
tational kinetic energy en = nr^Q^. 

We now use to represent in turn f, Q?, and en- We have 
then S = Si^if) and V5 = {dS j dM')l .^f . Rewriting Eq. Q 
in spherical coordinates (r, Q, (p) and knowing that in these coor- 
dinates the r- and 0-components of the effective gravity geff 
are 



or 

1 5<Dg 



r de 

we readily obtain 



+ Q''(r,6')rsin6'cos6', 



(10) 



— —[1-a^^+a^GM^] =0,(11) 

or r 06 

where we have used the notations 

1 dS 



and 



Cpd.ye 



(12) 



{ rsin^e if = / 
^ ^\ ifjf = Q2 

1 sin" 9 



(13) 
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The characteristic equation of Eq. (fTTT i is 
dr dd 



(14) 



As for the assumption S - S (J^) made here, and S are con- 
stant along the characteristic curves of Eq. [TT] and the differen- 
tial equation Eq. (fT4l i integrates immediately to give 



1 



sin^ 6 
r sin^ 9 
1 

sin^ 



= -Tj+Cj 
r 



if = 


/ 


if = 




if Jf" = 


en 



(15) 



with the C,^ as the integration constants (where has different 
meanings). Once the function giving the angular velocity Q.s{0) 
in the stellar surface is specified, the iso-rotation curves can be 
obtained easily everywhere inside the star using the relations in 
Eq. ( fTsT l. Each iso-rotation curve depends on an integration con- 
stant C,^ obtained as Cjsr = C,^(T^^, 0,,), where fl, is the colat- 
itude angle at which a given iso-rotation contour intersects the 
stellar surface. The constants Yj>r are defined as 



( 2a^GM \ 



(16) 



Since we consider convective regions where VS < 0, we find 
that Tjff > 0. We note that depending on the meaning of Jif, 
the expressions = Tj7?e, "^n/Rl and T^^ are dimensionless 
constants of the order of unity (Rg is the stellar radius at the 
equator). Defining 



1 (d\nPp-y 

y 



d\nr 



d\nr l\d\nJ^) 
which in convection regions is A,^ > 0, we have 



^ CM 



\i-A 



(17) 



(18) 



From th e model of a rapidly rotating star with M = 2QMq 
given m iMaeder et all (l2008l) . we derived the adiabatic gradi- 
ent V = dlnPp ''' /din r shown in Fig.[T] This quantity is shown 
as a function of stellar radius in the polar and equatorial direc- 
tions. The value of Y,^ depends also on the adopted value for the 
gradient dlnr/dlnQ? that represents the differential rotation. 
For the Sun, the term d In r/d In comes from seismology data 
and its value ranges from 1 to 10. However, until now in early- 
type stars, the factor li In r /li In has been a perfectly unknown 
quantity. For a rough estimate, we may assume that it can be at 
least of the same order of magnitude as in the Sun. Si milar values 
are als o obtained in the models calculated bv lMevne t & Maeder 
("200^. Using the values of the gradient V shown in Fig.[T] it can 
be shown that the ratio A^/{1 - A,^) is nearly constant over 
the entire convective region and that depending on the assumed 
value of dlnr/dlnQ?, it changes from to 1. Hence, we have 
that 



< Y,^ < 2(GM/rIq\ 



(19) 



where 0.5 < (R^Q} IGM) < 1 for rapid rotators. For the sake of 
comparison, we note that in the conv ective region of the Sun, it 
is |V| ^ 10 ^' and (GMIRl^l^) ^ 10^ (lBalbuJ2009l ). 



Examples of isorotation curves obtained from Eq. ( fTsT i are 
shown in Fig.|2l where we used stellar external contours whose 
rotational deformation is characterized by the ratio of centrifugal 
to gravitational acceleration at the equator rf,, - il^R^IGM - 0.8 
and a Maunder surface differential rotation law with parame- 
ter a - 0.3. The equation of the stellar surface with latitudi- 
nal differential rotation is discussed in some detail in Sect. [3] 
Solutions similar to those shown in Figs. and |2}3 were ob- 
tained by Balbus (2009) for the Sun. This author assumed that 
dPI dr » dS I do and dS 186 » dPI dr, which according to our 
expressions are for the limits a^f and anO? — > that suit 
slowly rotating objects like the Sun. 

2.3. Comments on the rotational profiles obtained 

The rotational profiles Q.{r, 9) in the convective layers beneath 
the surface of early-type rapidly rotating stars were inferred here 
using the marginal condition of the Solberg-H0iland stability cri- 
terion against axisymmetric perturbations. The derived solutions 
are the consequence of a dominant thermal wind balance where 
the entropy is: a) a function of the specific angular momentum, 
S -S(py, b) a function of the angular velocity, S -S(Q.^), and 
c) a function of the specific rotational kinetic ene rgy, S = S (en). 
The first case was widely studied in the literature dTassoulll 19781 
references therein), and a number of independent numerical at- 
tempts to obtain baroclinic stars inexorably ended up with re- 
sults close to the category shown in Figs.|2^1, i.e. conservative- 
like (lUrvu & Eriguchilll994 Il995h . Two-dimensional, implicit 
hydrodynamic numerical simulations by Deupree (1998) and 
iDeupred ('2001) showed that in convective regions O tends to 
adopt a cylindrical-like profile rather than a shellular or Q = 
constant form. Intermediate solutions to the classes shown in 
Figs.|2j3i and b2 are compatible with the seismological data of 
the Sun relative to its convective regions. The rotational profile 
obtained by Espinosa Lara & Rie utord (2007) for fully radiative 
stars curiously enters into the category of solutions shown in|2};2. 
We note that for all cases shown in Fig.|2]we have used constant 
values ofTjff all over the star. However, they are valid only for 
the convective regions, which lie roughly above the shaded sec- 
tor and where, as akeady noted, Y,^ is fairly uniform. 

According to a strong claim that shear instabilities widely 
prevail over all other mechanisms that can act to redistribute the 
angular velocity, most if not all models of stellar evolution of 
early-type stars with rotation were carried out by assuming that 
Q. has a shellular nature. This strongly simplifies the numeri- 
cal aspect of the study, but does not necessarily preclude that 
other distributions of O might coexist with the shellular one, 
or even dominate in wide domains of the stellar interior. The 
S = 5'(0^)-case with Yn > 1.6 clearly depicts shellular-like 
rotation over a wide interior stellar region. 

Although the present discussion should not be considered an 
argument in favor of non-conservative-like profiles of the an- 
gular velocity in the envelopes of early-type rapid rotators, the 
success that the displayed arguments do have in explaining the 
solar rotational profile in the convective layers, strongly moti- 
vate us to address the study of the surface differential rotation in 
rapidly rotating early-type stars. Whatever the observational or 
theoretical indication that dQ./d6 + can exist in the surface of 
rapidly rotating early-type stars, they would immediately imply 
that profiles of non-shellular rotation exist, perhaps of some type 
similar to that in Fig.|2] They would also immediately enable us 
to achieve clearer understandings of the transport of angular mo- 
mentum in the stars and the related mixing processes of chemical 
elements in the stellar surface layers. 
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X/Rp x/Rp 

Fig. 2. Q(r, 0)/Qe - constant curves inside the stellar envelope of a model star with a global shape generated with an equatorial 
rotational parameter rjo - 0.8 and a surface differential rotation represented by a Maunder relation with parameter a = +0.3. All 
curves were generated with parameters Tj ^.^jj - 0.3 and 1.6. The iso-rotation curves in (ai) and (ai) correspond to S - S{j^y, those 
in (bi) and (b2) are for S = S{Q?y, the curves in (ci) and (C2) obey the condition S - Sim^Q?). The solutions are supposed to be 
valid in the convective zones above the shaded domain. 



2.4. Induced gradient in Q(0) in the surface of early-type 
rapid rotators 



In rapidly rotating stars, there is a strong latitudinal temperature 
gradient, which can induce a latitudinal gradient on the surface 
angular velocity Qj. This effect can be estimated from Eq. (|4|l 
by calculating the latitudinal variation in the temperature over 
an isobar. Giving to the surface temperature a functional form 



T - T{P, 6), where P is the pressure and the colatitude, we 
obtain 



^ ' dP 86 r 



(20) 



where eg is the colatitude unit vector and r is the radial spher- 
ical coordinate. Making use of the equation of state P <x pT, 
Eq. | 2 0t and the r-component of Eq. (|5]l, Eq. (|4]i transforms into 
(lEspinosa Lara & Rieutordll2007h 



gr ( d\nT 
86 



(21) 
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where the notation (...)p indicates that the partial 0-derivative 
is along a barotropic surface. As in rapidly rotating stars, it is 
5 In r/50<O and gr<0, from Eq. (|2B it appears that dn/dz< 0, 
which means that D. should increase from the pole to the equator 
Although in rapid rotators there are convective zones in the enve- 
lope, the dominant radiation flux can still be estimated using the 
diffusion approximation Fiad = -X^T, where x is the coefficient 
of radiative co nductivity. F or simplicity, we adopt von Zeipel's 
approximation dvon Zeipel| [l924) to represent the stellar surface 
temperature distribution as a function of the colatitude 

T(e) = Tp{^j-^ , (22) 

where Tp is the polar temperature of the star, g = ^(gr)^ + (ge)^] ^ 
with gr and gg given by Eq. ( fTOb . and gj, is the value of the gravity 
in the pole. Using the Roche approximation for the gravitational 
potential, Eq. ( 1211 1 written in spherical coordinates becomes 



0.02 




^(radians) 

Fig. 3. Normalized surface angular velocity Q, as a function of 
the colatitude 8 and the rotational parameter 



dln^l^ tanfl^, ^dlnO? 

[l-(l-7?,.)/l(0,77,.)]- 



dr 



89 



--Me,71r)f2(0,r]r) 

r 

with the following definitions 



(23) 



I - T]r sin 



1 - T]r{2 - rir) sin 



2-77,.(l +rir)sv[?e 



sin 6* COS 6* 



c?lnr 



+ 2-/7, 



2 J 



GM 



(24) 



The characteristic equation describing the dependence of the 
surface angular velocity with 9 is then 



^ [ 1 - ( 1 - Tlr)fl (9, T]r)] \fx (9, rj,)f2{9, TJr) 



(25) 



which is solved by numerical integration. Figure [3] shows 
the 0-dependence of the angular velocity given as AQ(0) - 
[il(9,rir) - Q(40'',77,)]/Q(40'',77,-), which was obtained from 
Eq. dZSl) for different values of rj. In this calculation, the gradi- 
ent d In r/ 89 in the stellar surface was estimated using the Roche 
potential of a rigid rotator characterized by average rotational 
parameters 77. In Fig. [3] we see that in general it is 8Q./89 > 0, 
in particular for j],, > 0.7 near the equator. This acceleration 
seems to be similar in character t o the one previously calcu- 
lated bv lEspinosa Lara & Rieutordl(l2007l) . As 77 0, the curves 
AQ(0) ^ 0. 



3. Equation of the stellar surface 

3. 1 . The angular velocity distribution in the stellar surface 

Since the shape of a rapidly rotating star depends explicitly on 
its surface angular velocity distribution, we adopt a law to study 
the possible effects it can produce. In Sect. 12.41 we have seen that 
the temperature gradient induced by the gravity darkening effect 
introduces a small acceleration of the angular velocity from the 



pole towards the equator, which is not simple to represent analyt- 
ically. The physical properties of the layers beneath the surface 
could perhaps enforce a stronger surface angular velocity gradi- 
ent than inferred in Sect. 12.41 If in rapidly rotating massive stars 
rotational profiles of the class obtained in Sect. 12.21 actually ex- 
isted, a first approach to describing their surface rotation could 
rely on the use of a solar-like surface angular velocity. The so- 
lar surface differential rotation depends on the colatitude 9 as 
(Snodgra ss..l984) 



= 45 1.5- 65.3 cos^e- 66.7 cos'^e nHz. 

27r 



(26) 



This expression carries three coefficients that need to be deter- 
mined empirically. At the moment, it seems unrealistic to use a 
similar expression for other stars since their surface cannot be 
resolved. However, within the errors smaller than 1.6%, Eq. (l26b 
can be reduced to the simplest one 



^0(9) 
2n 



= 459.3(1 -0.29 cos^ 61) nHz, 



(27) 



where there are only two quantities to determine from observa- 
tions. A first inquiry about the differential rotation on the surface 
of massive and intermediate mass stars can then be justified by 
using the simplified expression (Maunder formula) 



Q(0) = Qo(l +acos^9) 



(28) 



In Eq. (l28l l. a < Q indicates that the pole rotates slower than 
the equator and vice-versa if a > 0. If massive and intermediate 
mass stars had surface rotations similar to the Sun, their differ- 
ential rotation parameter would be expected to be at least of the 

order of a 0.3. 

In general, to study stars with unresolved surfaces we should 
probably use expressions of the type 



Q(0) = a,[l+a/(0)], 



(29) 



where the function f{9) needs to be somehow specified in ad- 
vance. 
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Table 1. Radii ratios R^/Rp as a function of parameters ri„ and a 




Fig. 4. Shape of stars having average surface rotational parame- 
ters rjo = 0.7 (a) and 0.9 (b), and latitudinal differential rotations 
given by Eq. (l28T l for several values of the parameter a. The stel- 
lar shape for or = is indicated by a dashed line to more clearly 
show the effect of a in other cases. 



3.2. Formulation of the equation of ttie steiiar surface 

The shape of a rapidly rotating star is generally described us- 
ing the total potential of a self-gravitating system with a conser- 
vative rotation law, i.e., an internal cylindrical angular velocity, 
Q - Hivj), of which the rigid rotation is a particular case. When 
a star has a non-conservative internal rotational law O = Q(t(7, z), 
as used in models of stellar evolution and/or suggested by the 
developments given in Sect. 12.21 it is no longer possible to de- 
fine a rotational potential. However, according to Maeder (2001 
the surface of a star with a non-conservative rotation law can be 
identified as the region where an arbitrary displacement ds does 
not imply any work done by the effective gravity geff 



geff.liS = , 



(30) 



The effective gravity geff written in cylindrical coordinates has 
the form 



(31) 



where Og is the gravitational potential. Using the function 
*P(w, z) defined as 



(32) 



1o 


a = 0.6 


0.4 


0.2 


0.0 


-0.2 


-0.4 


-0.6 


0.4 


1.317 


1.275 


1.236 


1.200 


1.167 


1.138 


1.114 


0.5 


1.388 


1.339 


1.293 


1.250 


1.211 


1.176 


1.146 


0.6 


1.456 


1.401 


1.349 


1.300 


1.255 


1.215 


1.179 


0.7 


1.522 


1.462 


1.404 


1.350 


1.300 


1.255 


1.215 


0.8 


1.585 


1.520 


1.458 


1.400 


1.346 


1.296 


1.253 


0.9 


1.645 


1.577 


1.512 


1.450 


1.392 


1.340 


1.293 



the effective gravity can be expressed in the form 

geff = -[V^(C7, Z) + ^CT^VQ^] , (33) 

so that the condition given in Eq. ( l30b can be rewritten as 
d'i'im, z) + ^m\VQ.^.ds) = 0. (34) 

The function *F(ct, z) can be a potential if and only if it is 
dQ./dz - 0. It has been shown bv lMevnet & Maeded (11997 ) that 
in stars with shellular internal rotational profiles, the surfaces 
T(zrr, z) = constant are parallel to isobar surfaces. However, since 
d^'(m, z) is not an exact differential, the integration of Eq. ( l34l l 
depends on the chosen path. To define the shape of a star, it 
seems natural to integrate Eq. ( l34l i over a meridian curve, from 
the pole (spherical coordinate 6 = 0) towards an arbitrary point 
Rs(ff), where /?s(6') represents the function describing the stellar 
surface. Thus, by virtue of Eq. (|32] | we obtain 



^G{e)--n\m,zW + 



if 

2 Jo 



m^iVn^.ds) 



<I>g(0) . 



(35) 



which is in principle the sought equation to calculate the shape 
of a star having non-conservative rotational laws. Nevertheless, 
since the integration indicated in Eq. ( |35] | depends on the angu- 
lar velocity distribution defined only over the stellar surface, the 
following change in the coordinates describing the stellar surface 



mJ0) = R(e% sin ; z,(e) = R,(0) cos , 



(36) 



readily transforms VQ^(ziT, z).£/s into a total differential that en- 
ables us to integrate Eq. ( |35] | by parts and derive the sought re- 
lation 



1 , ldm^\ 
Og(.)--|q?(.)(— j^^.Oc 



(0). 



(37) 



In Eq. ( l37T i. we wrote to indicate explicitly that the angular 
velocity concerns only the stellar surface. From Eq. ( l37b . we can 
see that the stellar shape depends on both the rotational profile 
in the surface 0.^(0) and the gravitational potential Og, which 
in turn depends on the internal rotational law by means of the 
centrifugally distorted density distribution in the stellar interior. 
However, this last dependence is of second order or negligible, 
as will be shown in Sect. 15.41 



3.2.1 . Shape of stars with a Maunder rotation law in the 
surface 



According to the discussion in Sect. l2.4l and the simplified solar 
surface rotational law ( 1281 ). the differential rotation parameter in 
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fast rotating early-type stars is expected to be a < 0. However, 
because of the still poorly known physical characteristics of the 
convective external envelope layers of these objects, and pos- 
sible effects linked to magnetic fields, we cannot exclude that 
a > in some objects. In this respect, we note that hydrody- 
namic calculations in stars with 8.75M0 carried out bv lDeuprea 
disss') show that the rotation law in the convective core becomes 
of conservative type Q. ~ tir""^ (zcr being the distance from the 
rotation axis). If convection in the envelope forced it to rotate 
in a similar way, we should expect an external imprint of this 
law given by an a > 0. Many calculations of the stellar structure 
with fast differential rotation have been done in the past assum- 
ing that H0ilnad's criterion for dynamical stability is satisfied 
and that radiative viscosity has significant effects on the sur- 
face layers. The resulting rotation law implies that the specific 
angular momentum increases wit h the mass contained in cylin- 
ders with radius ht (lTassoullll978 ). which also leads to rotation 
laws that imply a > 0. However, the final shape of the rotation 
law in the envelope must certainly be the consequence of the in- 
terplay among many hydrodynamic and magneto-hydrodynamic 
instabilities, whose very final combined effect is still highly un- 
known. We adopt then the surface rotational law with the func- 
tional form given by Eq. dZTl ). and assume that the parameter a 
can imply either equator-ward or polar-ward accelerated angular 
velocities. 

In a first step, we separate the effects produced on the exter- 
nal stellar shape by the surface differential rotation from those 
induced by the internal rotation on the mass distribution. This 
can simply be done, as we justify in Sect. 15.41 by using a Roche 
approximation for the gravitational potential 



(t>c(0) = -GM/RM . 

Equation dJTl) can then be given the form 

RM ^ 1 

where we have written 



m 



if 

2 Jo 



If 



d0. 



(38) 



(39) 



(40) 



Since in relation ( [39l ) the function of the stellar surface Raiff) also 
appears in the integrand in terms of vr by ( |36] |. R^iff) is obtained 
by iteration. We start the iteration with a first estimate of linjl) 
based on the shape of a star with a rigid rotation characterized 
by the rotational parameter rj,, 



2d3 



^iR 



GM 



so that 



Rt\e) 













^ sin^ 6i| 



(41) 



(42) 



It is also obvious that in relation ( l39l l the equatorial-to-polar radii 
ratio Re/Rp is the solution of the iteration, which finally gives 



— = 1-1- rjoll —, — ,a 
Rp \ 2 



n R, 



(43) 



We recall that for rigid rotation in the Roche approximation we 
have 



RA , 1 

P / rigid 



(44) 



In Table [T] we present the ratios of equator-to-polar radii 
Re/Rp as a function of rjo and a used to obtain the stellar shapes 
shown in Fig.|4] which shows the shapes of stars computed with 
equatorial acceleration ratios rjo = 0.7 and 0.9, whose surface 
angular velocity is given by Eq. (l28T l. where the parameter a has 
several positive and negative values. In this figure, we can see 
that the stellar shape is a function of a that sensitively differs 
from that of homologous rigid rotators with the same parame- 
ter T],,. We note that for Q.{6) with a > (acceleration from the 
equator towards the pole), we always find that 



Rp(j]o, a) 



Rei^o) 



Rp(r]o) 



(45) 



rigid 



For modest a >0 ( recalling that in the Sun, it is a ^ —0.3 ), 
we can have Re/Rp > 1.5 even for rjo < 1, although the stel- 
lar surface has no polar dimples as in models with high rota- 
tional energy content (See Sect. I3.2.2I I. This is due both to a 
rotational stretching of and to the concomitant shrinkage of 
the polar radius Rp, where Qpoie > f^equator- When Qs(^') is ac- 
celerated from the pole towards the equator {a < 0), we have 
Reirjo, -\a\)/Rp{r]o, -\a\) < (/?e(?7o)/^p('7o))rigid- Here, the equa- 
tor/pole radii ratios respond mainly to the rotational stretching 

of Re. 

We note that in stars where a ^ 0, the Vsin / parameter does 
not only depend on the rotation in the equator. A straightfor- 
ward interpretation of the quantity Vsin / can produce an incor- 
rect value of the acceleration ratio rfo- Since we are dealing here 
with fast rotators {ri„ > 0.5), the formation of spectral lines used 
to determine the Vsin / needs to be treated properly by taking 
into account the surface velocity fields, the stellar deformation, 
and the gravity darkening effect, all dependent on t]o and a. We 
also note that if the interferometric measurements produce radii 
ratios Re/Rp > 1-5, it does not necessarily mean that the star is 
in a state of "supercritical" rotation. From these comments and 
the results shown in Table [T] we clearly see that the study of the 
surface differential rotation prefigures the need for two sources 
of information: a) spectroscopy, which can help us to estimate 
a and the inclination angle factor sin /; b) interferometry, which 
provides information related to the stellar geometrical deforma- 
tion and incidentally with a. We briefly discuss these issues in 
Sects. OandO 



3.2.2. Stars with polar dimples 

Owing to the centripetal acceleration, the fraction of the grav- 
itational force that must be supported by the pressure gradient 
increases towards the center of the stars. This means that the 
pressure gradient is more centrally concentrated in a rotating 
star than a non-rotating one. The rotation tends then to decrease 
the moment of inertia, which is accompanied by an increase in 
the density gradient towards the center. The gravitational attrac- 
tion towards the pole is thus stronger, which produces a flatten- 
ing of the pole or a dimple. In models of stars with high dif- 
ferential rotations in depth that enables the star to have a sig- 
nificant content of rotational energy, i.e. when the energy frac- 
tion T = {K = rotational kinetic energy; W = gravita- 
tional potential energy) becomes t > 0.03, the polar region is ei- 
ther strongly flattened or hollower dBodenheimed 1 97ll Clemeni 
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Fig. 5. Shape of model stars having the surface angular velocity 
profile given by relation ( l46l l for yS = 5.0 and several values of 
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Table 2. Models of stars with rigid rotation 



11974 Il979t IZoredll986l: iJackson et al.ll2005h . Similar dimples 
can also be produced by a surface angular velocity accelerated 
enough towards the polar latitudes. To illustrate this effect, we 
adopt the following angular velocity profile, which is sometimes 
used in the literature to m odel internal conservative (cylindrical) 
laws (iJackson et al ] l200l for stars with high rotational energies 
(0.02 < T < 0.22) 



1 +^[R,(0)/R^]^sm^e 



(46) 



where R^iff) is the equation of the stellar surface and y6 is a free 
parameter A relation similar to Eq. (l46T l may suit the rotational 
profiles shown in Fig. |2^ if the angular velocity is accelerated 
from the equator to the pole. Figure|5]shows the shape of model 
stars with a surface angular velocity profile given by expression 
(|46] | for j3 = 5.0 and several values of ratio t],, in the equator 
calculated with Eq. ( |39] |. It can be shown that using the Roche 
approximation and the rotational law in Eq. (l46T l. the ratio of 
equatorial to polar radii is then given by 



R, 
R. 



(47) 



Writing ^5(6*) = x + z , the condition for a flattened/hollowed 
polar region to appear is dz/dx > at jc = 0. The condition to 
obtain a hollowed polar region is then given by 



1 +rio 



I +13 



(48) 



Making u=Re/Rp = 1 +rio(l+/3)/2, relation (l48T l becomes w" ({u- 
1 ) > 2( 1 +y8). Since it is [m V(m - l)]min for m = 1 .5, condition (l48T l 
requires that /?> 2.375 independently of the value of rj,,. 

Since the simplest equation to represent the rotationally de- 
formed surfaces shown in Fig.|5]is 



RM 

Re. 



1 



1-^ 

Re 



COS 9 . 



(49) 



an alternative expression to relation (l28T l for the surface angular 
velocities accelerated from the equator to the pole could be 



Q,(0) 



1 +yS[l - (1 - |^)cos6i] ain^e 



(50) 
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1.297 


310 


4.377 


0.581 


0.081 


0.98 


9.265 


13.189 


1.369 


336 


4.522 


0.722 


0.087 


0.99 


9.266 


13.186 


1.410 


350 


4.571 


0.802 


0.088 


1.00 


9.267 


13.183 


1.511 


378 


4.619 


0.996 


0.090 


Note: Pl- 


is given in \ 


; cm ^ ; Q^-,. is given in s ' ; V^^. 


is given 


in km s ' ; the units of // M 


are cm" s 


; / is age and ^ms ifi 


the MS life span 









would require us to infer the ratio Rp/Re from interferometry, 
which according to Eq. ( l47T i is related to rjo and /3, and /3 derived 
from spectroscopy. However, before using expressions of the type 
given by either Eq. Ii46\l or (\50]l , it would perhaps be preferable 
to await observational confirmation that stellar shapes similar 
to those depicted in Fig. \5\actually exist. 

4. The gravity darkening effect in stars with 
non-conservative rotation laws 

Owing to the radiative equilibrium that prevails in the atmo- 
sphere of massive and intermediate-mass stars, the emerging 
bolometric radiative flux F is proportional to the temperature 
gradient Vr in the outer layers of the stellar envelope. Since 
these layers are also in hydrostatic equilibrium, their pressure 
gradient VP oc VT is balanced by the surface effective gravity 
geff. We have then the following phenomenological relation that 
describes the principle of von Zeipel's relation 



F oc Vr = VP oc geff 



(51) 
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As F cc T*^, where T^ff is the effective temperature at a given 
point 7 in the stellar surface, it follows the known von Zeipel's 
relation (von Zeipel 1924) 



,eff 



1/4 



(52) 



However, according to Poincare-Wavre's theorem dTassoull 

Il978l) only when the rotational law in the external layers is con- 
servative do the constant equipotential, density, and tempera- 
ture surfaces coincide to consider that "^(r) - - constant. 
Otherwise, as for rotational laws discussed in Sect. 12.21 ^€(T) is 
not constant with the colatitude Q. Its detailed expression for the 
non-conservati ve case with shellular rotation was obtained by 
lMaeded(fT99l . 

An equivalent relation to Eq. ( |52] | for spectroscopic and in- 
terferometric analysis can be given the form 



Xgeff' 



PcD 



(53) 



where Pgd - 0.25 + 6, so that 6 masks the variation in 
■^(r) over the stellar surface. The particular notation pcD 
used for the gravity darkening power in relation ( |53] | is in- 
troduced to avoid confusions with the fi parameter that ap- 
pears in the internal rotation law (l46T l. Two-dimensional (2D) 
models of rotating stars and some observations sugg est that 
< 5 < q [c-f. Eu cv (1967); Lovekinetal. (2006); van Bell eetaD 
(120061) ; iMomTier et all (120071) ; IZhao et al.1 (I2009D 1. We insist on 
6 not being in principle the function of taking into account the 
change in the stellar surface from radiative to possible convec- 
tive equilibrium-dominated layers because of the strong change 
in Teff with 6, as seems to be understood in Aufdenberg et al. 
(l2006l) and IZhao et aP (l2009b . but mainly the non-conservative 
nature of the rotational law in the external layers. 

5. Models of rotating stars 

To study the effect of the centrifugally distorted internal mass 
distribution on the gravitational potential in the surface of a star, 
to test the validity of the Roche approximation, it is enough to 
calculate models of stellar structure where only the global dy- 
namical aspects induced by the rotation are considered. To this 
end, we differentiate the primary effects produced by the rota- 
tion from those induced by the evolution. The primary thermo- 
dynamic effects due to the stellar evolution are taken into ac- 
count by making use of the relations between the pressure and 
density calculated with one-dimensional models of stellar evolu- 
tion without rotation. We assume thus that the changes produced 
on the barotropic relation P - P{p) at a given evolutionary stage 
of a star by the several instabilities and the diffusion of chemical 
elements unleashed through the stellar evolution by the rotation, 
have second order effects on the establishment of the dynamical 
equilibrium of the rotating star In principle, nothing prevents us 
from using the barotropic relations derived with models of stel- 
lar evolution with shellular rotation, but the results will probably 
not be more reliable. They all are calculated for low energy ra- 
tios T - KI\W\, which do not correspond to the model inputs 
regarding the rotational law and energy tested in the present at- 
tempt. 

Since we are not interested in the precise description of all 
non-linear time-dependent phenomena associated with the vis- 
cosity and internal flows in rotating stars, and because the total 
energy carried by the meridional circulation is small, our mod- 
els are axisymmetric, steady state, and circulation free . Owing 
to these assumptions and Poincare-Wavre's theorem dTassoull 



Il978h . our model stars behave as barotropes. We then adopt in- 
ternal rotational laws of conservative form, Q = Q(7ir), where tu 
is the distance to the rotation axis. According to the discussion in 
Sect l5.2l the conservative laws are expected to produce stronger 
stellar geometrical deformations for a given amount of rotational 
kinetic energy than the homologous non-conservative ones. 

For conservative rotation laws, the gravitational potential 
<l>(ziT, z) and the density distribution p{m, z) in the rotating star 
are simultaneous solutions to the hydrostatic equilibrium equa- 
tion 



and the Poisson equation 
AO = A-nGp, 



(54) 



(55) 



where (vj, 0, z) are the cylindrical coordinates with z containing 
the rotation axis, e^^ is the unit vector perpendicular to the z-axis, 
P is the pressure, and j - Q.m^ is the specific angular momen- 
tum. 

Equations d54b and d55] ) are solved with the adopted comple- 
mentary barotropic relation 



P^Pip) 



(56) 



The P = P(p) relations used in this work have a two-component 
polytropic character 



P = ap'>'" + bp''' 



(57) 



where the constants a, b, ja, and y/, were adjusted so as to: a) 
reproduce the pressure P^ and the density p^ in the center of the 
non-rotating star of given mass and evolutionary stage; b) ensure 
a continuous distribution of the pressure-density relation at the 
radius of the stellar core; c) obtain the correct ste llar mass at the 
stellar radius as tabulated bv jSchaller et al.l d 19921) for 1-D evolu- 
tionary models for the initial metallicity Z - 0.02. The function 
d57] i is continued in the stellar atmosphere by another pressure- 
density relation calculated by Castelli & Kurucz (2003) for stel- 
lar atmospheres as a function of the parameters (Teff, log^). 

The first order effects due to the stellar evolution are thus 
accounted for by the pressure-density relations at the center of 
the star and the dP/dp gradients. An additional term in relation 
d57] i could also in principle take into account the presence of the 
convective regions in the envelope induced by the rapid rotation, 
but this was not done in the present approach. The only rota- 
tional effect on the P = P{p) relatio n considered here is by means 
of the mass-compensation effect dSackmannlll970h . which in- 
creases the density pc at the center of the star. To this end, we 
iterated p^ until the nominal stellar mass M was attained. This 
iteration also implied that the central pressure f c was changed 
in accordance. 

The gravitational potential <l>(tir, z) is obtained by solving the 
Poi sson equation gi ven in Eq. dSST l with the cell-method adapted 
by IClementI (1 19741) to stellar structure calculations. The den- 
sity distribution pivr, z) is derived from the integrated form of 
Eq. ~ 



Jp, 



•p r'^ 

— = <I)(tir, z) - Oc + I D^(-nj)-nTdT!T , 
p, P Jo 



(58) 



where is the gravitational potential at the stellar center. Given 
a rotational law Q(tiT) and a barotropic relation in Eq. d56] l, the 
common solution to equations d54] i and dSST i is performed over 
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the entire space. The iteration of <1) and p is stopped when the 
highest density difference in the (cr, z)-space is max{6p/p) ^ 
10"*. In our iterations, the virial relation 6 - [2(K + U) - 
W]/\W\ - 0, where K = kinetic energy, U - internal energy, and 
W - total gravitational potential energy, is verified to better than 
5 «i 2x 10-4 in the ZAMS models, and (5 * 6x lO^^ by the TAMS 
models. Since in the framework of conservative rotational laws, 
the surfaces of constant pressure, density, and total potential are 
parallel, the rotationally distorted shape of our models is defined 
by the total equipotential surface that contains the polar 'photo- 
spheric' radius Rp. This radius is identified by the layer whose 
density satisfies the model-atmosphe re relation tros s(p) - 2/3 in 
the stellar atmosphere models of Castelli & Kuruczl (.2003). The 
local effective temperature in the pole needs to satisfy also the 
gravity darkening ef fect. We mo dified accordingly the effective 
temperature given bv lSchaller et al. ( 1992) for the g iven mass M 
using von Zeipel's approximation dvon Zei pel 1924). The trans- 
formation to the rotation-dependent effect ive temperature wa s 
performed following the procedure given in lFremat et al.l(l2005h . 



5.1. Rigid rotation 

The simplest conservative rotational law is that of rigid rotation. 
The models of stars with rigid rotation were calculated assuming 
that at each evolutionary phase they experience instantaneous to- 
tal redistribution of their internal angular momentum. The char- 
acteristics of these models are given in Table |2] In this table, 
f/f MS is the fractional age of the star, where Tms is the time that 
a non-rotating star of mass M spends on the main sequence, Qcr 
is the critical angular velocity, Q/Qcr represents the angular ve- 
locity for which the model was calculated, p^ is the core den- 
sity of the rotating object, Rs/Rq is the equatorial radius of the 
model star in solar units and Rs/Rp the equatorial-to-polar radii 
ratio, Veq is the equatorial linear velocity in kms"', J/M is the 
total specific angular momentum, rj - Q.^RI/GM is the ratio of 
centrifugal to the gravitational acceleration in the equator, and 
is the ratio of the kinetic rotational energy (K) to the ab- 
solute value of the gravitational potential energy (W). Some of 
these models are shown in Fig. |6] where in all cases the iso- 
density surfaces are for the same p/p^ density ratios. Similar 
models were also calculated for other masses and age ratios 
///ms- In spite of the simple approach used to calculate them, 
the radius ratios co mpare very closely with those calculated by 
lEkstrometai] (1200 8) for the same energy ratios. We note, 

however, that the models obtained by these authors cannot be 
compared directly with ours because the distribution of the in- 
ternal rotational velocity is not the same. Starting from a quasi- 
rigid rotation in the ZAMS. l Ekstrom et al.l (12008 ) accounted for 
a consistent evolution of the angular momentum distribution in- 
side the star throughout the calculated stellar evolution span. 



5.2. Sheiiuiar differential rotation 



IZahn et aP (l2010h studied the shapes of stars with internal shei- 
iuiar differential rotation and concluded that in 7Mq stars the 
radius ratio Re/Rp at critical equatorial rotation can be enlarged 
from 1 % to 4%, depending on the internal rotational energy con- 
tent and the evolutionary phase. However, it can easily be seen 
that when two rotation laws are described with the same func- 
tion and both imply the same central to surface angular velocity 
ratios Qcenter/fisurface, but onc is shcllular and the other cylindri- 
cal, the law that is shellular may have a weaker effect on the 
global internal density distribution than the cylindrical one. For 




Fig. 8. Schematic comparison of the centrifugal force produced 
by a shellular and a cyUndrical rotation law when both are de- 
scribed by the same analytical function 



simplicity, we change the independent variable P (pressure) in 
the shellular law by the radius r. The comparison is depicted in 
Fig. [8] the specific centrifugal force acting at a point p in the 
stellar interior is 



^shellular 



^cylindncal ^ q2(^^ g)^ 



(59) 



for shellular and cylindrical angular velocity distribution laws, 
respectively (m = rsinO is the distance to the rotation 
axis). At the equator, if we assume that Q'=5''''"*"'='''(zi7,7r/2) = 
Q,sheiiuiar(^^ ;r/2), for a given point p(z, vj), it will be 



cylindrical 



^cylindrical^^^ ^ ^shellular ^^-j 



(60) 



This holds in particular for the points in the stellar surface. 
However, since it has been shown in Sect. 13.21 that the exter- 
nal geometrical deformation of the star is mainly a function 
of the angular velocity law on the surface, cylindrical rotation 
laws will carry stronger geometrical deformations than shellular 
ones for a similar total kinetic energy. Moreover, in Sect. l2.2l we 
have anticipated that strong cylindrical components of the an- 
gular velocity distributions in the external layers probably exist. 
It is then important to determine the order of magnitude of the 
quadrupole factor y in rotating stars using conservative internal 
rotational laws. To this end, in the next section we calculate two- 
dimensional models of rotating stars with conservative internal 
rotational laws. 

5.3. Cylindrical differential rotation 

The models with differential rotation were calculated using the 
same barotropic relations P-P{p) as for the solid-body rotation. 
We adopted the internal angular velocity distributed as 



Q(nT) = 



1 +/3{m/R,r ' 



(61) 
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Fig. 6. Rigid rotators, a: Iso-density surfaces in model stars of M - 5Mq in the ZAMS rotating at (ai) Q = 0; (a2) Q = 0.85Qc; (a.3) 

Q = Qc = 1.92 X lO""* s"'; b: Iso-density surfaces in stars of M = 5Mq in the TAMS rotating at (bi) Q = 0; (bi) Q - 0.85Qc; (bs) 
Q. - Q.C - 4.79 X 10"^ s '.The iso-density surfaces are labeled with the corresponding density ratios p/pc, which are the same in all 
panels of the figure. 




10 15 



10 20 30 40 50 60 70 

x/R,, 



10 20 30 40 50 60 70 80 

x/R,,, 



Fig. 7. Differential rotators , a: Iso-density surfaces in model stars of M = 15Mq in the ZAMS having the internal rotation law (1611 1. 
whose parameters Qg/^cr and p are indicated. Q.cr is for the critical rigid rotation. The iso-density surfaces are labeled with the 
corresponding density ratios p/pc, which are the same as in Fig. |6] and for all panels of this figure, b: Iso-density surfaces in model 
stars of M = ISM© in the TAMS having the internal rotation law ( |6TI ), whose parameters Qg/^cr and /3 are indicated. In each panel, 
the ordinates are in the same scale as the abscissas, but they differ from one to the other panel. 



13 



Zorec et al.: Spectroscopy and interferometry of early-type differential rapid rotators 



Table 3. Models of stars with internal differential rotation 



ZAMS 






M = 5Mo 




n„ = 1.92X lO-* 






ZAMS 






M = 15Mo 




n„ = 1.34X 10-" 






n„/n„ 




R,JRo 


RJRo 








J/M 


K/\W\ 




Pc 


R.v/Ro 


RJRo 




Ra/Rp 




J/M 


K/\W\ 
















xlO-" 


















xur" 










P = 


2 
















= 


2 










1.0 


19.996 


2.671 


2.71 1 




1.052 


117 


3.326 


0.0044 


1.0 


5.945 


4.915 


4.986 




1.047 


151 


9.920 


0.0064 


2.0 


21.593 


2.821 


2.967 




1.269 


225 


6.723 


0.0175 


2.0 


6.482 


5.092 


5.383 




1.265 


294 


20.091 


0.0250 


3.0 


24.741 


3.451 


3.938 




1.947 


290 


10.437 


0.0389 


3.0 


7.580 


5.828 


6.666 




1.878 


396 


31.115 


0.0553 


3.5 


27.451 


4.813 


5.799 




3.143 


255 


12.898 


0.0532 


3.5 


8.529 


7.330 


8.938 




2.809 


381 


38.020 


0.0752 


3.9 


32.121 


9.933 


11.652 




6.998 


152 


17.775 


0.0683 


4.0 


10.930 


16.655 


20.092 




7.455 


216 


58.322 


0.1035 


1.0 


19.870 


2.639 


13 = 
2.656 


4 


1.023 


70 


2.824 


0.0033 


1.0 


5.899 


4.875 


fi = 
4.915 


4 


1.025 


91 


8.220 


0.0046 


2.0 


21.022 


2.667 


2.729 




1.138 


139 


5.642 


0.0130 


2.0 


6.273 


4.861 


4.981 




1.133 


181 


16.443 


0.0179 


3.0 


23.095 


2.731 


2.867 




1.353 


202 


8.458 


0.0285 


4.0 


8.085 


4.955 


5.455 




1.697 


341 


32.810 


0.0663 


4.0 


26.378 


2.906 


3.185 




1.728 


250 


11.318 


0.0489 


5.5 


11.219 


5.774 


6.887 




2.843 


396 


45.878 


0.1163 


5.0 


31.538 


3.400 


3.899 




2.489 


269 


14.450 


0.0736 


6.2 


14.319 


7.599 


9.291 




4.558 


349 


56.026 


0.1457 


6.2 


46.905 


6.692 


7.790 




6.456 


181 


22.434 


0.1128 


6.4 


16.610 


9.221 


1 1 .247 




5.958 


304 


64.434 


0.1588 


1.0 


19.796 


2.634 


J3 = 
2.652 


6 


1.022 


50 


2.482 


0.0027 


1.0 


5.874 


4.856 


4.889 


6 


1.010 


65 


7.103 


0.0035 


2.0 


21.022 


2.626 


2.670 




1.095 


100 


4.946 


0.0104 


3.0 


6.681 


4.726 


4.869 




1.221 


196 


21.149 


0.0303 


4.0 


23.095 


2.631 


2.768 




1.426 


195 


9.761 


0.0392 


5.0 


8.684 


4.544 


4.941 




1.693 


324 


34.451 


0.0770 


5.0 


28.193 


2.677 


2.907 




1.727 


235 


12.103 


0.0586 


7.0 


12.871 


4.618 


5.417 




2.664 


424 


46.715 


0.1336 


7.0 


40.196 


3.160 


3.641 




2 969 


275 


17.048 


0.1047 


8.0 


16.731 


5.171 


6.299 




3.689 


430 


53.793 


0.1648 


8.35 


67.403 


5.514 


6.446 




6.938 


196 


25.828 


0.1483 


8.4 


19.300 


5.767 


7.106 




4.564 


408 


58.318 


0.1796 


1.0 


19.746 


2.626 


/s = 

2.640 


8 


1.009 


39 


2.227 


0.0022 


1.0 


5.857 


4.857 


/3 = 
4.884 


8 


1.011 


51 


6.293 


0.0029 


3.0 


21.808 


2.580 


2.630 




1.164 


118 


6.610 


0.0190 


2.0 


6.092 


4.790 


4.837 




1.066 


102 


12.548 


0.0112 


5.0 


26.464 


2.519 


2.669 




1.495 


195 


10.776 


0.0492 


4.0 


7.143 


4.551 


4.722 




1.304 


209 


24.704 


0.0425 


7.0 


34.953 


2.529 


2.788 




2.070 


264 


14.665 


0.0878 


6.0 


9.316 


4.229 


4.610 




1.733 


319 


35.796 


0.0869 


9.0 


50.57B 


2.873 


3.343 




3.305 


292 


18.741 


0. 13 13 


8.0 


13.338 


3.980 


4.573 




2.373 


428 


45.301 


0. 1363 


10.0 


66.836 


3.572 


4.219 




4.962 


264 


22.199 


0.1575 


9.6 


18.901 


4.078 


4.940 




3.327 


484 


52.380 


0.1764 


TAMS 






M = 5Mo 




n„ = 4 


79 X 10-* 






TAMS 






M = 15Mo 




n„ = 2.76 X 10-" 






1.0 


27.626 


6.639 


/3 = 
7.708 


2 


1.051 


71 


1.692 


0.0008 


1.0 


9.255 


13.564 


= 
13.777 


2 


1.049 


86 


4.1 19 


0.0007 


2.0 


28.085 


7.556 


7.915 




1.280 


132 


3.482 


0.0031 


2.0 


9.357 


15.546 


16.370 




1.285 


160 


8.521 


0.0030 


3.0 


29.015 


12.077 


13.648 




2.334 


137 


5.633 


0.0072 


3.0 


9.572 


25.858 


25.259 




2.411 


162 


14.021 


0.0073 


3.5 


29.954 


27.310 


30.655 




5.467 


78 


7.475 


0.0104 


3.5 


9.982 


66.541 


73.661 




6.333 


81 


19.758 


0.0107 


1.0 


27.601 


6.535 


J3 = 
6.564 


4 


1.028 


43 


1.541 


0.0006 


1.0 


9.249 


13.363 


= 
13.483 


4 


1.026 


52 


3.753 


0.0006 


3.0 


28.632 


7.844 


8.169 




1.384 


110 


4.829 


0.0059 


2.0 


9.331 


14.200 


14.476 




1.128 


99 


7.644 


0.0025 


4.0 


29.723 


10.200 


1 1.025 




1.984 


116 


6.805 


0.0107 


3.0 


9.480 


16.212 


16.970 




1.387 


133 


11.882 


0.0058 


5.0 


31.729 


19.386 


20.867 




4.110 


81 


9.748 


0.0177 


4.0 


9.732 


21.574 


23.300 




2.021 


138 


16.979 


0.0108 


5.5 


34.006 


42.148 


44.030 




9.199 


43 


13.124 


0.0221 


5.0 


10.236 


44.869 


48.918 




6.602 


87 


25.481 


0.0183 


5.7 


36.217 


85.056 


88.533 




19.108 


22 


16.649 


0.0246 


5.5 


11.028 


129.78 


136.36 




9.820 


35 


39.665 


0.0240 


1.0 


27.584 


6.497 


13 = 
6.515 


6 


1.021 


31 


1.428 


0.0006 


1.0 


9.246 


13.246 


fi = 
13.307 


6 


1.013 


37 


3.479 


0.0006 


3.0 


28.440 


7.200 


7.383 




1.232 


83 


4.396 


0.0051 


3.0 


9.437 


14.815 


15.268 




1.238 


72 


10.791 


0.0050 


4.0 


29.275 


8.086 


8.447 




1.494 


100 


6.025 


0.0091 


4.0 


9.626 


16.812 


17.571 




1.498 


121 


14.905 


0.0091 


5.0 


30.516 


9.898 


10.573 




2.007 


103 


7.883 


0.0145 


5.0 


9.914 


20.975 


22.325 




2.043 


123 


19.763 


0.0148 


6.0 


32.471 


14.334 


15.413 




3.190 


88 


10.297 


0.0216 


6.0 


10.393 


31.844 


34.264 




3.419 


99 


26.539 


0.0223 


7.0 


36.522 


32.588 


33.300 




7.737 


48 


14.903 


0.0308 


7.0 


11.634 


91.051 


95.665 




10.725 


42 


43.227 


0.0334 


1.0 


27.571 


6.474 


f3 = 
6.496 


8 


1.018 


24 


1.337 


0.0005 


1.0 


9.243 


13.239 


J3 = 
13.327 


8 


1.006 


29 


3.260 


0.0005 


3.0 


28.311 


6.933 


7.037 




1.173 


67 


4.082 


0.0045 


3.0 


9.408 


14.238 


14.472 




1.164 


82 


10.008 


0.0049 


5.0 


29.992 


8.280 


8.640 




1.610 


94 


7.099 


0.0127 


5.0 


9.790 


17.301 


18.070 




1.624 


113 


17.656 


0.0128 


6.0 


31.368 


9.791 


10.361 




2.081 


96 


8.858 


0.0185 


6.0 


10.112 


20.816 


22.169 




2.148 


113 


22.314 


0.0189 


8.0 


36.726 


20.408 


21.171 




5.150 


65 


14.213 


0.0345 


7.0 


10.602 


28.159 


30.082 




3.159 


99 


28.363 


0.0267 


9.0 


47.921 


88.304 


88.136 




24.993 


18 


24.175 


0.0453 


8.5 


12.570 


86.558 


90.257 




11.410 


41 


49.916 


0.0438 



Note: Pe is given in g cm 
nally deformed object 



; £lcr is given in s ' ; Veq is given in kms ^ ; the units of J/M are cm^s ^ ; R^v/Re is the radius of a sphere having the same volume as the rotatio- 



where /? > is a free parameter and is the angular veloc- 
ity at the axis of rotation. For simplicity, we considered n-2. 
Nevertheless, it can be shown that all models with n < 2 obey the 
Solberg-H0iland stability criterion, and that at whatever distance 
m the energy ratio Tim) is K{'nT)l\'W{m)\ < 0.138. This ensures 
that no region in the stellar interior is unstable to the secular in- 
stability that can carry it to a three-axial Jacoby ellipsoidal con- 
figuration. 

In Table [3] we indicate the characteristics of the calculated 
models. The entries have the same meaning as in Table |2l The 
higher values given in Table [3] are about the highest 

we could obtain for the given masses, ages, and jS parameters. 
Figure |7] shows some of these calculated models. 



Using the same algorithm, we can obtain models that ap- 
proach more closely the limit of dynamical stability ~ 1/4 
if the objects are on the ZAMS and have masses higher than 
M - 15Mq. However, Clement ( 1979) in his analysis of the sec- 
ular stability of differentially rotating objects already pointed out 
the unstable character of model stars with > 0.1. 



5.4. The Roche approximation 

With the models thus calculated, which give us the shape of a 
rotationally distorted star, we can test to what degree the gravi- 
tational potential in the stellar surface may deviate from the 
simple central-field expression used in the Roche approximation. 
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Table 4. Quadrupole factors in stars with rigid rotation at differ- 
ent angular velocity ratios Q/Q^ 





n/f2, = 


0.5 


0.7 


0.9 


0.95 


1.0 




Phase 






r 






5.0 


ZAMS 


0.0012 


0.0018 


0.0026 


0.0028 


0.0030 


5.0 


TAMS 


0.0003 


0.0007 


0.0016 


0.0020 


0.0025 


15.0 


ZAMS 


0.0015 


0.0025 


0.0037 


0.0040 


0.0043 


15.0 


TAMS 


0.0004 


0.0007 


0.0018 


0.0023 


0.0028 



iHubbard et al.l il97 ^ showed that the gravitational potential 
of rotating centrally condensed objects can be given by a multi- 
pol expansion 



<l)r 



GM 

"m 



n=I 



(62) 



where R,, is the radius of the rotationally undistorted star, 72« are 
the zonal harmonic coefficients, and f 2«(cos 6) are th e Lege ndre 
polynomials. To derive relation ( |62] |. iHubbard et al.l (1 19751) as- 
sumed that the density can be expanded in Legendre polynomi- 
als. The J2,, coefficients are then obtained as integrals over the 
volume of the axisymmetric object which generally produce val- 
ues so that J2 > 0, J4 < 0, 76 > 0, 78 < 0, and so on. Since we 
are not interested in developing a detailed theory of the gravita- 
tional potential of rotating stars, but only to test the validity of 
the Roche approximation, we do not calculate the volume inte- 
grals over the rotationally deformed star, but use the harmonic 
coefficients coefficients as mere parameters giving a quantita- 
tive indication of the deviation from the central-field form of 
the gravitational potential. We determine the quantities 72b in 
expression ( |62] | by simply fitting the gravitation potential pre- 
dicted by the 2D models of rotating stars calculated in Sects. ISTI 
and l5.3l using a least squares method. In this algebraic estimate 
of 72;,, their order of magnitude are preserved, but not always the 
expected sign. 

In relation ( |62] |. we retain only the quadrupole moment to 
estimate the effects in rigid rotators 



1>r 



_gm( 

m\ 



Ro 



R(0) 



72P2(C0S 0) 



(63) 



The shape of the surface thus becomes 

-.2 



Re 



(64) 



where we have introduced the notation 



7 - Ji 



(65) 



The largest difference between the estimate given by Eq. ( |64] | 
to the ratio Rsi6)/Re and that produced by the Roche approxima- 
tion in Eq. (l44l i. is expected for = 0. Thus according to ( l64l i we 
have 



R^ _ 1 + j{r]o + y) 
Rn ~ 



(66) 



Since it is always true that Rg/Rp > 1, relation ( 1661 ) implies that 
for whatever < rjo < I and y ^ the equator to polar radii 
ratio is slightly larger than obtained from the sheer Roche ap- 
proximation. 

In TablelH we indicate the quadrupole factors y for 5Mq and 
15Mo stars in the ZAMS and TAMS evolutionary phases ob- 
tained by fitting relation ( l64l l with surface gravitational potential 
obtained with the model calculation. In this table, we see that the 
smaller 72 the more evolved is the star, simply because the star is 
more centrally condense d such that t he central-field approxima- 
tion for holds better. IZahn et alj (12010) calculated the same 
factors for rigid and shellular differential rotators. Reducing their 
72 estimates by (Rp/Rg)^ to be able to approach the ratio (Ro/Re)^ 
more closely than using the average stellar radius for Rq, we 
see that for the ZAMS and TAMS epochs of rigid critical ro- 
tator s we obtain the sam e values for the harmonic coefficient 72 
than lZahnetaLl (1201 Oh . Having then y < 0.004, from dUli we 
see that neglecting the quadrupole term in (l64l l. in rigid rota- 
tors at critical rotation the Re / Rp can be underes timated by less 
than 2%. For shellular rotators. ' Zahn et al] (1201 0) found that at 
critical equatorial rotation the flattening changes by 4% because 
^^center/f^surface = 4 in the ZAMS and is smaller near the TAMS. 

Finally, we note that within approximation ( l64b . the actual 
ratio of the centrifugal to gravity acceleration at the equator 
should be 



riJ = 



i + ir 



(67) 



In stars with conservative differential rotation laws of the type 
given in Eq. dFTT l that imply rotational energies t = K/\W\ > Tc,-, 
the deformations of the internal mass distribution and the devi- 
ations to the central-field gravitational potential are expected to 
be stronger than in stars with rigid rotation. To test the use of the 
Roche approximation in these cases, we calculated the zonal har- 
monic coefficients y2n = JiniRo/Re)^ from « = 1 to 4. Table|5]lists 
the zonal multipolar coefficients y2„ obtained by fitting relation 
(|62] | with the calculated gravitational potential for several val- 
ues of the coefficient /? and angular velocity ratio Qo/Qcr, where 
Q.0 is the angular velocity at the rotation axis, and Q„ is the 
critical angular velocity of an homologous star (same mass and 
evolutionary stage) behaving as a rigid rotator. We insist on the 
algebraic nature of the y2n coefficients, since they were derived 
as fitting parameters to model-calculated gravitational potential. 
In this Table, we also give the fractional deviations of the radius 
vectors AR(0)/R(0) = [/?(0)™dei - R(0)yJ/R(e)^oA.i x 10^ for 
different colatitude angles 0, where R{0)modei is obtained with 
the models calculated in Sect. 15.31 while the R(0)y2„ values re- 
fer to those estimated with the total potential that includes the 
gravitational potential given by Eq. i62\ . A quick inspection of 
Table |5] shows that in most cases a fairly precise description of 
the stellar surface can be obtained by considering only 72 and 74 
in the gravitational potential given by Eq. ( |62] |. Nevertheless, in 
all cases the Roche approximation can be valid to better than 5% 
for ratios ilcenter/^^surface - 4. The approximation works better 
the more evolved the star is and the higher the rotational energy 
parameter t. 

6. Attainable stellar parameters 

According to Vogt-Russell's theorem, "the complete structure of 
a star in hydrostatic and thermal equilibrium is uniquely deter- 
mined by the total mass M and the run of chemical composition 
throughout the star, prov ided the structure e quations are func- 
tion of local parameters" dCox & Giulilll968h . Then, for a given 
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chemical composition, denoted here in a generic way by 2, it 
follows that parameters that conventionally specify the physical 
properties of a star are the luminosity L, the mass M, and the 
radius R (or the stellar age t). 

An object in permanent rotation can still be considered to 
be in hydrostatic and thermal equilibrium. However, owing to 
the numerous effects the rotation induces, the number of param- 
eters needed to characterize an object can be high: mass, total 
luminosity, age, total angular momentum, internal distribution 
of the angular momentum, dimensions or stellar geometry, in- 
clination angle, degree of non-uniform distribution of the tem- 
perature in the surface, mixing of chemical elements, and degree 
of their non-uniform internal and surface distribution, etc. This 
number is certainly higher than we can tackle and/or determine 
from observations. We consider here only parameters we can ac- 
tually determine with the diagnostic tools we presently have at 
our disposal. However, since many of these quantities are model- 
dependent, the resulting characterization of rotating stars would 
probably be limited or incomplete. 

We may classify the accessible stellar parameters into three 
classes: 1) those related to the stellar internal structure; 2) pa- 
rameters relevant to the structure of the atmosphere over the ob- 
served stellar hemisphere; 3) quantities related to the apparent 
stellar geometry. 

Models of stellar structure : The existing models of stars with 
rotation start their MS evolution with low rotational energies in 
the ZAMS, i.e. 



T = KI\W\ < t':«(M) , 



(68) 



where K is the rotational kinetic energy, W is the gravitational 
potential energy, and r"/ is the energy ratio for a rigid rotator at 
a critical equatorial rotation in the ZAMS. For stars with masses 
ranging from 3Mq to some 6OM0, it is on average t"/ < 0.015 
(IZor ec et al. 1988a). Fom Table |2] we see that tzams is a func- 
tion of M and Veq. For rigid rotators in the ZAMS, Veq and Q/Qc 
are equivalent parameters. 

Owing to mass and angular momentum losses, the energy ra- 
tio T becomes a function of the time, so that t(M, f) < t'^%^^s{M). 
On account of possible redistribution phenomena of angular mo- 
mentum that could take place in the pre-MS evolutionary phases, 
we should characterize the stellar initial rotational law with at 
least one single quantity. This single parameter should also cor- 
respond to the expected value of the kinetic energy ratio when 
Tzams > "'"zams ci - However, existing models of rotating stars 
that provide us today with evolutionary paths to infer the stellar 
masses and ages, are calculated with initial rigid rotations in the 
ZAMS. In this case, t, Q/Oc, and Vgq can be assumed to carry 
the same information. We must keep in mind that the masses 
and ages we infer from these models are depend on the assump- 
tion given in Eq. ( |68] |. The models provide us with qu antities 
averaged over the ste llar surface (cf. iMevnet & MaedeJ (i2000l) ; 
lEkstrom et al.l (l2008l) l namely surface-averaged bolometric lu- 
minosity (L), effective temperature (Teff ), surface effective grav- 
ity (loggeif) (occasionally), and equatorial linear rotational ve- 
locity Veq 



(L) = {L{M,t,Q, tzams)}, 

{T,ff} = {T,niM,t,Q,TzAMs)) , 

<log geif ) = (log gsdM, t, Q, Tzams )) , 

Veq = Veq(M,f,e, Tzams) ■ 



From iMaeder & Mevnetl (l2000l) . we can see that because 
of rotation, stars behave as though they have lower effective 
mass than in reality, which mainly affects the emitted bolo- 
metric luminosity . According to this mass-compensation effect 
(ISackmannI 1 1 970h and the core fuelling with hydrogen by the 
meridional circulation, rotating stars have longer MS life-spans. 
Owing to the change in the opacity of the envelope by the dif- 
fusion of chemical elements from the core, mainly He, stars be- 
come more luminous than their mass-homologous rotationless 
objects. Evolutionary timescales in the post-MS phases are also 
sensitively changed. If stars were actually neat differential ro- 
tators, Le^_wdthrote^ energies tzams > T"r (iBodenheirneil 
ll97ll : IClementl[T979 ). the evolutionary tracks needed to interpo- 
late stellar masses would also differ from those existing today. 

Models of stellar atmospheres : High resolution spectra and 
spectral energy distributions of rapidly rotating stars are both 
/-dependent (/ is the inclination angle). They can be modeled 
by taking into account the rotationally-induced stellar geometri- 
cal deformation, but mainly the concomitant latitude-dependent 
surface gravity and temperature d istribu tion: gravity darkening 
(lFrematetalJl2005l: iLovekin et al. l l200 6': 'Gillic h et al.ll2008h . 

Since a priori we do not have any indication on what may be 
the function describing the surface angular velocity of a star, in 
Sect. l3.1l we adopted a simple Maunder-like expression given by 
Eq.|28] i.e., Q(0) - Q(Qo, a, 6), where Q„ is the angular veloc- 
ity of the equator, a is the surface differential rotation parame- 
ter, and is the colatitude angle. Instead of using O,, we prefer 
to use the force ratio parameter rj,, defined in Eq. (HD) . Since 
a ?i 0, the iso-radial velocity lines responsible for the Doppler 
shifts and causing the rotationally broadening of spectral lines 
are no longer straight lines but curves, so that the true velocity 
in V sin / does not necessarily corre spond to the equatorial linear 
velocity dZorec et al.lll988bll2004l) . To this adds the uncertainty 
in the V sin / parameter related to the less effective contribution 
to the line broadening b y the eq uatorial regions affected by the 
gravitational darke ning (fstoecklevi 1 968bl: pTownsend et al.l2004l ; 
iFremat et al.ll2005h . 

For the colatitude-dependent distribution of the effective 
temperature, we concluded in Sect.|4]that Tgff{0) -^i Y. geft'^™, 
where can be assumed to be caused by the average stellar geo- 
metrical deformation and where pcD - 0.25 +d has a free param- 
eter 6 that should b e ^ < according to what is suggested in the 
literatur e [cf. Lucvi (1 1967h; ILovekin et al. (2006); van Belle et aTl 
(I2006h ; iMonnier et al.1 (l2007h ; IZhao etlll (l2009l) 1. For models 
that begin evolving in the ZAMS as rigid rotators, t and r\o are 
also synonymous. The apparent, hemisphere-dependent bolo- 
metric luminosity, effective temperature, effective gravity and 
rotation parameter are then functions of at least the seven un- 
knowns of mass M, age t, initial chemical composition Q, equa- 
torial rotation parameter 770, surface differential parameter a, 
power pGD in the von Zeipel relation, and the incUnation angle / 



L"PP = L(M, f,e, 770, a, pGD,0, 

r7 = reff(M,f,e,/7„,a,|3GD,0, 

log = log ^(M, f, e, ^o, a, Pgd, , 

(VsinO'P'' = y(M, f, e, T/o, a, Pgd) X sin / 



(70) 



(69) 



We also include the dependence on Pgd in log^'^PP and 
(Vsin /)"PP, since their prediction is based on spectral lines pro- 
duced by gravity-darkened models. The above seven-fold de- 
pendence is valid also for the equivalent widths (W,\) of spec- 
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tral Hnes, the line FWHM (full width at half maximum), the ze- 
roes q„ of the Fourier transforms of spectral lines, as well as the 
energy distributions and spectrophotometric parameters such as 
{Ai,D) of the BCD system, which describe the Bal mer discon- 
tinuit y and are highly sensitive to (logg^PP, T^JJ'') dZorec et"al] 
l2009l) . Referring to all these observable quantities with a generic 
letter ff^, they are formally 

= ^,(M,f,e,77„,a,PGD,0. (71) 

We note that for stars with a given initial chemical compo- 
sition Q and a prescription for their internal rotation, the para- 
metric couples ((L), (reif )), and (M, t) are equivalent quantities. 
From observations, we derive L^pp, T^^, and log g^PP, which are 
/-dependent, while the models of stellar evolution provide us 
with (L), (Teff), and sometimes (log^). This means that M and 
/ can be estimated by iteratio n, where we also iterate simultane- 
ously the inclination angle / (IZorec et al. 2005). 

Stellar geometry : Since it is possible to use the Roche ap- 
proximation to describe the rotationally induced stellar deforma- 
tion, from Eq. dJTl l it comes that the inclination angle-dependent 
shape of rapidly rotating star depends only on rj,,, a and /. If we 
are not interested in carrying out imaging of fast rotating objects, 
interferometry can provide us with data that to a first approxima- 
tion do not depend on M, t, Q, and |3gd, i e. the apparent shape 
of the star (y, x), where x is the coordinate measured along the 
equatorial line projected onto the sky, and y is: 

R(rio, a, /, x) 

y{x) = , (72) 

SO that the apparent polar-to-equatorial radius ratios are y(x - 
0) - Rp{i)/Re and y(x = 1) = 0.0. We examine this in greater 
detail in Sect. 16.21 

If models are calculated for a given metallicity Z, which 
stands for the initial generic quantity Q used above, we conclude 
that a massive and intermediate-mass fast rotating star with sur- 
face differential rotation can be described as a function of six 
parameters 

M, rjo, t, a, Pgd, and / , (73) 

which we should be able to derive by interpreting the observed 
data. We briefly review some of the data and techniques that can 
help us to determine some of these parameters, in particular j]o, 
a, and /. 

6.1. Information on the differential rotation from spectral lines 

As stated in Sect. |4] the distribution of the effective temperature 
in the stellar surface deviates from von Zeipel's prescription, i.e. 
Tetf oc^jj™ with pGD ^ 1/4 if Of 7^0. However, in what follows of 
the present discussion we shall assume the validity of the origi- 
nal von Zeipel relation, but use a modified expression of the ef- 
fective gravity geff according to the latitudinal dependence of the 
surface differential angular velocity 0.(7], a, 0) given by Eq. ( l28l l. 

In this section, we report first estimates of the effects pro- 
duced by surface rotational laws Oij], a, 9) on some spectral lines 
and envision with them the way we can obtain estimates of a and 
the inclination angle / from spectroscopy . 

We modified our FASTROT code (iFremat et al.ll2005h . ini- 
tially written for a uniform angular velocity on the stellar sur- 
face, to take into account an external Maunder-like rotational 



profile and the related changes in the gravity and the effective 
temperature as a function of the colatitude in the calculation 
of spectral line profiles. In Fig.|9^, we show the Sin 6347 line 
profiles obtained with the new FASTROT program, which are 
produced by MS stars with T^s - 19500 K and \ogg - 4.0 at 
rest, whose rotation is characterized by several ratios of the an- 
gular velocity in the equator Oo/Qc without surface differential 
rotation (a - 0). In Fig.|9}3, we show the Sin 6347 line profiles 
for stars with the same {T^s, log g) rest parameters as before, but 
now for Qo/Qc = 0.9 and several degrees of differential rotation 
{a + 0). All line profiles in Figs. |9^ and|9]3 were obtained with 
gravity-darkened models seen at the inclination angle is / = 45". 
We note that the kinematic effect on the line broadening with 
surface differential rotation is determined by curves of constant 
Doppler shift, which are not straight lines as in the solid body 
rotation, but curves who se shapes depend on the value o f a and 
the inclination angle / (Dom iciano de Souza et a l. 2004b). Other 
glimpses into the sensitivity of different spectral lines in early- 
type rapid rotators to the effects carried by a surface rotation law 
tl(r\, a, 0) can be seen in lZorec et al.l (l2007h . 

The Fourier transform can give a quantitative description of 
the effects induced on the spectral lines by the latitudinal differ- 
ential rotation and perhaps on the associated gravity darkening 
effect. This information can be carried by the ratio qilqx of the 
second (^2) to the first zero {q\) of the Fourier lobes, as sug- 
gested by Reiners & Schmitt (2004, and references therein). In 
Fig.fTOh. we see the ratio ^2/^1 determined by the Fourier analy- 
sis of the lines shown in Fig.|9j), as a function of the ratio Qo/Qc 
for stars without surface differential rotation {a - 0). We see 
in this figure that the ratio qilqi becomes sensitive to the grav- 
ity darkening only when Og/Oc approaches a critical value. In 
Fig.fTOb. we show the ratios ^2/^1 for the profiles in Fig.|9j5, i.e. 
Qo/Qc = 0.9 and a Q, where the corresponding gravity dark- 
ening is also taken into account. As in Fig. |9] the inclination in 
Fig.[Tn]is / = 45". From the results shown in Fig.fTOb. we see that 
the differential rotation in fast rotating early-type star produces 
fairly large effects on the ratio qilqi, which are then expected to 
be measurable. It is also expected that a series of relations of the 
type 

92/^1 =/(«,/) (74) 

can be obtained, whose functional dependence on a and / differ 
from line to line. This is because the dominant excitation mech- 
anism of spectral lines is not the same for all of them, and that 
they can react more or less in different ways, or be more or less 
sensitive to the non-uniformity of the surface effective tempera- 
ture and gravity. This precludes our being able to define a proto- 
type line-broadening function that could be used for all spectral 
lines analyzed in a spectrum, which will be studied in detail in a 
following paper. 

We note that for this type of study, only spectral lines that do 
not undergo strong Stark broadening can be used. Accordingly, 
in early-type stars the Hei 4471 and Mgii 4481 lines are more 
freque n tly used to deter mine th e Vsin ; parameter (c .f. lSlettebakl 
(ll982h : IChauville et al.l (l200lh : lFremat eTaP (l2005h '). To test the 
sensitivity of these lines to the effects related to the gravita- 
tional darkening, we calculated the line profiles of Hei 4471, 
Mgii 4481, and Sin 6371, and obtained the diagrams shown 
in Fig. [TTh . where for each studied line the inclination of the 
rotation axis is presented against the equivalent width and the 
FWHM (full width at half maximum). For the calculations 
shown in Fig. [TT] we set the differential rotation parameter to 
a-Q. However, similar behaviors occur when ai^Q. The curves 
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Fig. 9. Rotationally broadened Sill 6347 model line profile in a main sequence B2-type star (rest Teff = 19000K, log § = 4.0) seen 
at the inclination angle i-45°. (a): Line profiles affected by gravitational darkening in stars rotating at different surface O/Qc rates 
and differential rotation parameter 0.0; (b): Line profiles with gravitational darkening effect in stars rotating at Qo/Qc = 0.9 and 
for several values of the differential rotation parameter a. To calculate the line profiles, we used (3gd = 0.25. 



in Fig. [TTIshow the differentiated sensitivity of the Hei, Mgii, 
and Si ii lines to the effects carried by the rotation. 

Using both the measured equivalent widths and FWHM of 
lines, and the diagrams of the type shown in Fig. [TTh . we obtain 
intersections as shown in Fig. [TTb that determine estimates of 
Q/Qc and From Eq. (l74l i. we can then infer a value of a and 
use it to obtain a new series of curves similar to those in Fig. [TTh. 
This iteration can be followed until convergence for Q./Q.C, i, and 
a is attained. To infer a first estimate of uncertainties related to 
this method, we note that the curves of Fig.[TT]can be reasonably 
fitted with relations such as 



(75) 



where oj - Q/D.^. It then follows that uncertainties 6W/W in the 
equivalent widths carry errors in the inclination angle / which 
are on the order of: 



Si 

- - Bico) 
I 



1 - 



C(w) 



6W 



(76) 



For an order of magnitude estimate, we consider Achernar, for 
which measurements of the apparent radii ratio Re/Rn were ob- 
tained (Domician o de Souza et alJ 

l2003tlKervella & Domiciano de Souzai200 6) and another analy- 
sis of the inclination angle was performed ( Vinicius et al. 2006). 
Assuming for this star that rjo ~ 0.8 - 0.9 and / < 75", from 
the Hei4471 and Mgii4481 we can infer, respectively, that 
B4471 = 3.3 and C4471 = 11.6, B4481 = 2.5 and C4481 = 13.7, and 
6W/W ~ 0.03 that lead to 6i ~ 5", which is an upper limit, since 
as seen in Fig. [TTb. the Si 11 lines reduce this uncertainty. 

A simple polynomial fit of a as a function of the ratio 
= li/li shown in Fig. [TO]leads to uncertainties of 6a ^ 0.15 
for Rq = 1.65 and <5/?q - 0.1, and smaller values of 5a for 
other ratios 7?q. However, the determination of the parameters 
Q/Qc (or rjo) and /, and then a, can be greatly improved us- 
ing curves similar to Fig. [TT]calculated for many other spectral 
lines of different excitation character, i.e, collision-dominated or 



photoionization-dominated, which enab le then to infer the incli- 
nation angle factor sin / more precisely (iHutchings & StoecklevI 
1977; Ruusalepp 1989). A future paper will be devoted to the 
discussion of this method and the related uncertainties. 

6.2. Information on the differential rotation from 
interferometry 

Interferometry can help provide different types of information on 
rotating stars: a) data related to the geometrical shape of stars; 
b) independent estimates of the differential rotation parameter a 
and the inclination /; c) imaging of the non-uniform hemisphere 
of a rotationally deformed star. In this introductory paper, we do 
not address the imaging techniques. 

a) Stellar geometry. The main goal is to measure the ge- 
ometrical deformations carried by the stellar rotation, and 
in particular those added by the surface differential rota- 
tion described in Sect. 13.2.11 Using the integrated fluxes 
in different spectral bands, the geometrical shape of the 
sources can be determined from the visibility measure- 
ments using the wi dest possible coverage of the plane 



(u,v) ( van Belle et al.l [ 2001; Domiciano de Souza e t alj 2003 



Aufdenberg et al. 2006; van Belle et al. 2006; McAHster et al 
2005; Kervella & Domiciano de Souza 2006; Monnier et al 
2007; Zhao et al. 2009). The resulting description of the geome- 
try of the source is however i-inclination angle dependent. 

As in Eq. ( [72l i. we write as y(x, i) the coordinates measured 
perpendicular to the x axis drawn in the direction of the stellar 
equator projected onto the sky. The y(x, i) ordinates determined 
from the interferometry are then related to the true shape of the 
star by 



y(i,x) = [Rjie, r]o, a) sin^e-x^]^'^ cos i + 
Rsi0, rjo, a) cos sin / , 



(77) 



where Rf^iO, rj,,, a) is the equation of the stellar surface normal- 
ized to /?e given by Eq. ( [39l l and is the colatitude angle mea- 
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Fig. 11. (a): Inclination angle as a function of the equivalent half-intensity width of the Hel 4471 line for gravitationally darkened 
stars rotating at the indicated Q.ol0.c values, where the differential rotation parameter is set to a = 0.0. ; (b): Intersection of constant 
FWHM(Hel 4471)- and FWHM(MgII 4481)- curves of an early-type rapid rotator seen nearly pole-on. 



sured in the stellar reference system. Equations (fTTT i and (|39T l 
are two independent equations that for a given pair of prameters 
{i, a) enable us to obtain an x = x{ri, ff) relation and determine 
the sought-after parameter rj. This leads us to the final objec- 
tive, i.e., to derive the true radii ratio Rp/Re = R^iO - 0)/Re of 
the star from the apparent stellar shape y{i,x)/Rf.. To achieve 
this purpose, the inclination angle / can be inferred using several 
iteration and modeling methods: 1) approximate ones such as 
those tested by Vinicius et al. (2006) and Carciofi et al. (2008); 
2) the spectroscopic inferences described in Sect. 16.11 3) spec- 



trophotometric constraint s imposed by the ( Ai , D) parameters of 
the Balmer discontinui tv (IZorec et al.ll2005h: 4) far and near-UV 



spectrophotometry ( Fremat et al .l l2002Ll2005h : 5) differential in- 
terferometry as described below (see item b). 

A first insight into the uncertainties incurred when using the 
VINCI/VLTl interferometric data can be inferred by assuming 
that the apparent star corresponds to the project ion of an ax- 
isymmetric ellipsoid. Equation dTTl l then becomes (iKanaan et al.l 
[2008 ) 



Pap = [l -(1 -pl)sin^i] 



1/2 



(78) 



where we have 

Pap = (-Rp//?e)apparent and ptr = (RJR.),n,.. FoT the VINCWLTI 



data on Achernar analyzed by iKervella & Domiciano de Souzal 



(l2b0 6). we infer that Sp^p - 0.02. Using a Monte Carlo simula- 
tion for the propagation of errors, we obtain Spa - 0.03, which 
from Table [1] implies that the differential rotation parameter a 
can be estimated with an uncertainty of 6a<Q.15. 

Improved estimates o f Pap are expec ted from the 
VEGA/CHARA instrument dMourard et al .' 2009). To determine 
the lowest stellar fattenings that we can reliably measure and the 
related uncertainties, the important quantity we need to estimate 
is the visibility ratio 



y2 _ y2 

' p ' e 



(79) 



where is the visibility in the polar direction of the apparent 
shape o f the star and is the visibility in its equatorial di- 
rection. IMourard et alJ (l2009l) demonstrated that highly reliable 



measurements with the VEGA/CHARA instrument can be per- 
formed of the angular diameter when the relative uncertainty in 
the squared visibility curve is crv2/V^ < 0.02. We then expect 
that visibility curves obtained in directions towards the appar- 
ent stellar polar and equatorial directions, respectively, be able 
to be differentiated and thus a difference between R^ and Rp be 
found as soon as dV^/{V^} > 0.02. We then tested the ratio in 
Eq. ( 1791 ) for BO, B2, and B8-type objects of luminosity classes 
V and III with apparent magnitudes V = 5 mag and V - A 
mag. The magnitude y = 5 mag corresponds to an upper con- 
servative limiting magnitude for the VEGA/CHARA interferom- 
eter in the medium-resolution spectral mode. The tested B sub- 
spectral types were chosen to be BO to include massive objects 
for which rapid rotators are frequent, B2 spectral type which rep- 
resent the highest frequency of Be stars, and B8 to represent Bn 
stars, which are rapid rotators without emission lines and most 
currently found among late B-type stars. 



Us i ng the fundamental pa ramet ers given in IZorec & BriotI 
(119911), iDivan & Zore^ (|1982|) . and Zorec et al. (2009) for the 
visible absolute magnitude My, the absolute bolometric mag- 
nitude Mboi, and the effective temperature T^s, respectively, 
we inferred the circular angular diameters 6o of the tested B 
sub-spectral types associated with their apparent average radius 
Ro - (/?e+^p)/2. We calculated the visibilities for two baselines, 
B = 156 m and B = 300 m, attainable with the VEGA/CHARA 
instrument at the effective wavelength A - 5500 A, assuming 
that stellar discs have uniform brightness. The visibilities were 
calculated separately for the polar and equatorial directions and 
test radius ratios ranging from Re/Rp - 1.05 to R^/Rp = 1.25, 
to account for the lower deformations of rotators characterized 
by parameters (770, a) reported in Table [1] The results obtained 
are displayed in Table |6] We note that the fundamental parame- 
ters given in Tableware not for specific stars, but correspond to 
hypothetical objects having the quoted MK spectral types. 

From the results displayed in Table |6] we conclude that pro- 
vided the objects are observed with large enough baselines, even 
a modest flattening Re/Rp ~ 1.05 could be detected at the 
limiting magnitude V = 5 mag. Brighter objects with slightly 
larger flattenings are found to have visibility ratios dV^/{V^) 
that are easily beyond the uncertainty range of 3cry2 /V^ ~ 0.06. 
Hence, rotational parameters {T]o,a) can be inferred with the 
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Fig. 10. Ratio q2/q\ of the second to the first zero of the Fourier 
transform of the spectral lines shown in Fig. |9] (a) Ratios for 
line profiles that correspond to stars rotating at different Qo/Qc, 
which are affected by gravitational darkening, but without sur- 
face differential rotation (a = 0.0). (b) Ratios affected by differ- 
ential rotation and gravitational darkening, as a function of the 
differential rotation parameter a in stars where Qo/Qc = 0.9. All 
relations shown in this figure are for the inclination angle / = 45°. 



VEGA/CHARA instrument to much better than 2% in rotators 
brighter than V = 5 mag with rates ri„ > 0.4. 

For the sake of completeness, we derive in Appendix lAl an 
expression for the visibility ratio dV^/{V^) as a function of the 
apparent flattening Re/Rp and the argument x = 7:6o(B/A), where 
6o = (2/D)[(/?e + ^p)/2] is the average angular diameter, given 
by 



(V-) ^ \ \+Sl,[nx)IV{v)f 

T{x) = dV{x)ld\nx 

"ep - (RJRp)+l 



(80) 



which gives a clear insight into the dependence of dV^/{V^) on 
Rs/Rp, provided that and Vp are both in the respective first 
lobes of Ji(x), i.e. Xp < 3.8327 and Xg < 3.8327. In Fig. [12] we 
compare the ratios dV^/{V^} given as a function of the apparent 
flattening Rs/Rp and calculated for several average arguments x 
(circles). In the same figure, we also superimpose the estimates 
derived from the approximation given in Eq. ( l80t (dashed lines). 
Figure[T2]may help us to sharpen the observation strategy. From 
the figure, it is clear that when R^ ~ Rp, it is dV^/{V^) —> and 
when for the same baseline one of the visibilities is significantly 
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Fig. 12. Calculated (circles) and analytically approximated 
(dashed line) visibility ratios dV^/{V^} as a function of the ap- 
parent radii ratio R^/Rp, for different parameters x = nOoiB/A). 



higher than the other the visibility ratio saturates, dV^/{V^} —> 2, 
so that the resolution for the determination of high flattenings 
becomes less obvious. 

Finally, we note that for rjo - 0.4 and a ~ the flatten- 
ing is Re/Rp - 1.2, which also implies that QoOc - 0.63 and 
Veq/Vcri t - 0-7. This transla tes into Veq ~ 26Q[{M/Mo)/Ro/Ro)] 
km s^' (iFremat et al.ll2005h . which means that the equatorial ve- 
locities range from 300 to 220 km s"' in dwarf stars with spec- 
tral types going from BO to AO, and from 230 to 180 km s"' 
in gia nt stars of the same spectral typ e. Searching in the cata- 
log bv lGlebocki & Stawikowskil(l2000l) for emission-less objects 
hotter than A2, with Vsin; > 200 km s apparent magnitudes 
V<5 mag, and declinations 6 > -20°, we found 63 OB stars that 
can be observed with VEGA/CHARA in the framework of the 
present study. Since their average inclination angle is probably 
/ > 52° ~ arcsin(7r/4), we may expect that many of them have 
true equatorial velocities in the range 250 km s"' < Veq < Vcrit- 

A detailed account of the iteration procedure used to deter- 
mine Re/Rp, T]o, a, and / from interferometric data will be given 
elsewhere for actual stars observed with VEGA/CHARA. 

b) Parameters {a, i). IChelli & Petrovl 

(1995a''b^ and'Petr ov et all (Il996in iighlighted the advantages of 
combining Fouri er transforms with the differential interferom- 
etry. iDomiciano de Souza et alJ (l2004bh showed that using the 
wavelength-dependent photometric barycenter, measured in an- 
gular units of interferometric intensity maps, it is possible to 
produce independent observational estimates of the differential 
rotation parameter a and the inclination angle /. The method is 
fairly independent of the characteristics of the observed spectral 
line and gives well defined indications of the sign of a. It also 
provides two estimates of /. Combining then data from the above 
class a) with either interferometric or spectroscopic information, 
it i s possible to infer the right value of /. The method developed 
by IDomiciano de Souza et akl (j2004bj) is used for slow rotators. 
Its implementation to gravity-darkened and or-dependent differ- 
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Table 6. Visibility ratios dV'^HV'^) calculated for different B 
sub-spectral types and luminosity classes, with several equator- 
to-polar radius ratios Re/Rp. The visibility ratios are calculated 
atA- 5500 A for two baselines and apparent visible magnitudes 
V. 







B 


= 156 m 


B = 300 m 
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A 

45 / 
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u. lyj 

0.123 


U.U30 
0.015 


\J.\J\JJ 

0.026 


0.131 
0.052 


0.135 
0.050 


0.094 
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1 lA 
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U.ZDZ 

0.165 


U.UDj 
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0.044 


0.242 
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0.277 
0.096 


0.504 
0.176 
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0.362 










B2111 
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479 


0.250 
0.158 


0.056 
0.021 


0.107 
0.041 


0.218 
0.084 


0.248 
0.085 


0.453 
0.161 


0.868 
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57 
91 


0.433 
0.273 


0.036 
0.015 


0.067 
0.028 


0.139 
0.054 


0.145 
0.055 


0.269 
0.099 


0.536 
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115 
182 


0.414 
0.261 


0.173 
0.062 


0.315 
0.116 


0.627 
0.240 


1.329 
0.274 


1.750 
0.499 


1.966 
0.946 



Note: D is the distance and is the angular diameter. For a given 
spectral type, the first line is for V = 4 mag and the second for V = 5 
mag. 



ential rapidly rotating early-type stars is currently being devel- 
oped. 

c) Imaging. Recent imaging of fast rotators using interfero- 
metric data has detected deviations from the standard von Zeipel 
grav ity-darkening prescription (Monnier et al. 2007; Zha o et aP 
12009 ). The imaging technique can therefore provide comple- 
mentary information to solve a possible parameter degeneracy 
that may include the differential rotation in the stellar surface. 

We note that while in the imaging techniques parameters 
such as pGD and the inclination angle / are inferred mostly from 
confidence levels in the^^ fitting surfaces, in the present attempt 
we try to infer a and / using direct methods. This does not pre- 
clude, however, the use of confidence levels to decide the choice 
of other parameters, in particular the gravity darkening exponent 

pGD- 

6.3. Method to derive the parameters 

We developed GIRFIT, a power ful tool to fit obse rved spec- 
tra of low and high dispersion (iFremat et al.l 12006*) with the- 
oretical fluxes based on model atmospheres calculated with 
TLUS TY (Hubenv & Lanz 1995) and/or with ATLAS9 (Kuruc3 
Il993t ICastelli et al.. ,1997.) . Adopting parameterized values of 
Pgd, we can thus infer the apparent quantities given in Eq. ( TTOb . 
Corrections for first order rotation effects can then be attempted 
as in Frem at et al, (2005), Zorec et al. (2005), Martavan et al. 
( l2006h . and lMartavan et alj (l2007b to derive the average param- 
eters in Eq. ( |69] | and obtain estimates of M and f. Since in most 



remaining calculations the mass enters as M ' , we can proceed 
to derive the other quantities with approximate values of M. 

In principle, from Eq. (l74l i and relations of the type shown in 
Fig. [TTb we can derive a) and /. By minimizing the differences 
between Eq. (iTTl l and interferometric observations, another set 
of independent estimates of (z/o, a, i) can be derived. 

The choice of the iterated parametric set (M, f, rjo, a, i) as a 
function of the gravity darkening power Pgd can be finally de- 
cided by applying a GIRFIT-like residual minimization proce- 
dure between the observed and the FASTROT model spectra. 

7. The rotational profile in the envelope 

From the discussion above, it follows that if c 9^ 0, it will proba- 
bly be impossible to know anything more about the characteris- 
tics of the rotation law in the convective layers beneath the stel- 
lar surface. Apart from the rigid or shellular-type rotation pro- 
files being able to be excluded, spectroscopic and interferometric 
data do not provide any further information. Higher radii ratios 
Re/Rp than allowed by the inferred values of a may be indica- 
tive of large values of the energy ratio t - K/\W\, but no detailed 
description can be derived about the internal rotation law itself. 
To establish the properties of this law near the surface, at least 
whether they are to be likely of the speculated types generated by 
the condition S = S(j^), S = S (O^), or 5 = 5 (en), no information 
can be obtained otherwise than from the analysis of non-radial 
pulsation modes excited in these regions, provided they can be 
unamb iguous l y det ected from time series of spectral line vari- 
ations. lAndol (119801) showed that the "effective depths" of p„- 
modes, where n = 1 - 5, for stars with masses from 5Mq to 2OM0 
range from 0.7 to 0.97?e, which means that the detected pulsa- 
tions, mainly the splitting of modes, may carry some information 
on the characteristics of the rotation law in the upper envelope 
layers. In all cases, the interpretation of these modes rely on the 
knowledge of the stellar fundamental para meters, in particular 
the inclination angle ; of the rotation axis (Floauet et al. 20001 
[2002; Neiner et al J 120021: iLevenhagen et al. 2003; Fremat et all 
2006j) whose determination can be attempted with the methods 
outlined in the present paper 

8. Conclusions 

Recent theoretical works suggest that rapidly rotating early-type 
stars should have rather deep convective layers in the envelope. 
This situation may favor some coupling between the convective 
region and the differential rotation, as happens in the Sun. We 
have considered the possible properties the iso-rotation curves 
can have in the envelope immediately under the stellar surface 
by assuming several conditions for the specific entropy 5 as a 
function of the: 1) specific angular momentum j, as implied by 
marginal Solberg-H0iland stability condition; 2) squared angular 
velocity O^, which accounts for the Solar differential rotation in 
the convection zone; 3) specific kinetic energy ■nP'Q?, which re- 
produces the calculated rotational profiles in radiative envelopes. 
In all cases, the function that describes the gradient dQ./d6 in the 
stellar surface must be specified in advance. We thus assumed 
a simple Maunder formula with a unique differential rotation- 
parameter a to represent the imprint on the stellar surface pro- 
duced by the differential rotation beneath the surface. In spite 
of the simplicity of this relation, for many theoretical reasons it 
would be highly interesting to infer from observations reliable 
orders of magnitude of a and its sign. While a value a ^ can 
imply that the rotation beneath the surface is neither rigid nor 
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shellular, neither spectroscopy nor interferometry can say any- 
thing about their actual functional nature beneath the surface. 
Other piece of information is thus needed, perhaps non-radial 
pulsation data of pulsation modes excited in the intermediate and 
upper stellar layers. 

In early-type fast rotators, we have shown that the surface 
temperature gradient might induce a positive gradient in the sur- 
face angular velocity, i.e., dQ/dO > 0. 

In this paper, we have summarized the kind of information 
that it is possible to extract from the spectroscopic and inter- 
ferometric data about the differential rotation of rapidly rotating 
early-type stars. The equation of the surface of a star with a non- 
conservative rotation law was discussed. It has been shown that 
a differential rotation in the surface induces measurable stellar 
deformations. These deformations carry a surface gravitational 
darkening effect that needs to be studied consistently with the 
induced geometrical deformations. We have shown that the ef- 
fects induced by the surface differential rotation can be studied 
consistently and in a reliable way by a combined analysis of 
spectroscopic and interferometric data. From spectroscopy, es- 
timates of the differential rotation parameter a with uncertain- 
ties 6a < 0.15 and of the inclination angle ; with 6i < 5" can 
be obtained. The interferometry can provide information about 
the true ratio Re/Rp of the equatorial to the polar radii to within 
less than 3%, provided that a reliable estimate of the inclination 
angle is given from spectroscopy. The VEGA/CHARA interfer- 
ometric instrument is capable of providing useful information 
about rotationally induced stellar deformations for nearly 60 OB 
rapidly rotating stars (Vsin ; > 200 km s ' ) whose spectra are not 
marred by circumstellar emission/absorption. Differential inter- 
ferometry can help distinguish the a and i parameters and pro- 
vide reliable independent estimates of each. 

By calculating simpUfied models of stars at different evo- 
lutionary stages in the MS with internal rigid and differential 
conservative rotation laws, we have tested the use of the Roche 
approximation to represent the gravitational potential of rotating 
stars. We have concluded that we can safely use the Roche ap- 
proximation for most, if not all stellar objects. Only in highly 
centrifugally deformed objects may this approximation not be 
applicable, but doubts exist as to whether these cases actually 
exist in Nature. 

Uncertainties and constraints carried by the spectroscopic 
and interferometric modeling as well as the application of these 
methods to real stars will be presented in separate papers. 
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Appendix A: Approximate analytic representation 

of dV^KV^) 

We assume that the observed visibility of an oblate source having 
an equator to polar radii ratio Re/Rp to be given for an equivalent 
uniform disc with angular diameter 6o as 



V(x) 



2Ji(x) 



(A.l) 



where Ji{x) is the Bessel function of the first kind, x - jt6o{B/A) 
where B is the baseline of the interferometer projected onto the 
sky and A is the effective wavelength at which observations are 
carried out. To achieve an accuracy of better than 0.2%, we can 
use the development 



y(x)- 1 + ^, 



(A.2) 



n=l 



with ai = 0.0173261, 02 = -0.160464, and 03 = 0.229252, 
which is valid for 0.0 < x < 3.8327 [first lobe of Ji{x)]. Defining 
the argument x as a function of the average angular diameter of 
the source 60 = 2x[{Re+Rp)/2]/D (where D is the distance of the 
source), the corresponding arguments of the visibility function 
in the equatorial and polar direction of the same source can be 
written as 



Xe = X X (1 - Sep) 

(RJRp) - 1 



(RJRp) + 1 ■ 



(A.3) 



Since 0.0 < 6sp < I we can use a convergent development of 
Eq. (IA.2b in a series of powers of dep. Conserving only the first 
order terms in d^p, Eq. ( IA.2b becomes 



V^six) - V{x) ± 4pr(x), 



(A.4) 



where r(x) - 2^=i na„x" - dV(x)/d\nx. Writing Vp - V[x x 
(1 + ^ep] and Ve = V[x x (1 - d^p], we obtain 



dV^ 



' p ' e 



<V2) i(yp2 + Vi) [ 1 + 6lp[r(x)/V(x)] 



4p[r(x)/y(x)] 



(A.5) 



which is assumed to be used only when both and Vp are in the 
respective first lobes of 7i(x), i.e. Xp < 3.8327 and Xg < 3.8327. 



23 



Zorec et al.: Spectroscopy and interferometry of early-type differential rapid rotators 
Table 5. Zonal harmonic coefficients 72/! and residuals ARi6)/R(ff) 
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33 


40 


46 


49 


51 




8.0 


0.167 X H)-' 


-0.722 X 10-3 


0.296 X 10-" 


0.623 X 10-" 


-1477 


-1329 


-450 


130 


54 


23 


34 


41 


45 


47 




8.4 


0.155 X 10-' 


-0.620 X 10-3 


0.205 X 10-" 


-0.319 X 10-" 


-3411 


-3002 


-1074 


303 


137 


7 


30 


40 


45 


47 
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Table 5. continued. 







72 


74 


76 


ys 


0° 


5° 


15" 


30° 


40° 


50° 


60° 


70° 


80° 


90° 


8 


1.0 


-0.144 X 10"^ 


-0.726 X 10-^ 


0.267 X 10-" 


-0.402 X 10-' 


11 


11 


13 


20 


31 


39 


43 


45 


47 


49 




2.0 


0.107 X 10-^ 


-0.598 X 10-3 


0.167 X 10-' 


-0.300 X 10 ' 


7 


8 


12 


21 


32 


40 


43 


45 


48 


49 




4.0 


0.922 X 10"^ 


-0.464 X 10-' 


0.261 X 10-' 


-0.413 X 10-" 


-20 


-17 


4 


22 


33 


40 


44 


47 


50 


51 




6.0 


0.177 X 10"' 


-0.790 X 10-' 


0.237 X 10-' 


-0.358 X 10-" 


-81 


-72 


-23 


21 


36 


41 


46 


50 


53 


54 




8.0 


0.229 X 10"' 


-0.110 X 10-2 


0.781 X 10-" 


-0.450 X 10-5 


-243 


-219 


-64 


24 


33 


40 


46 


53 


57 


58 




9.6 


0.219 X 10-1 


-0.106 X 10-2 


0.488 X 10-" 


-0.119 X 10-5 


-1380 


-1237 


-394 


119 


56 


32 


43 


51 


56 


58 












M = 15Mo (TAMS) 




















2 


1.0 


0.833 X 10-^ 


0.230 X 10-2 


-0.373 X 10-' 


0.111 X 10-2 


-33 


-33 


-44 


-73 


-113 


-144 


-158 


-167 


-176 


-182 




2.0 


0.797 X 10"^ 


0.105 X 10-2 


-0.102 X 10-2 


0.344 X 10 


39 


32 


-11 


-71 


-116 


-1412 


-157 


-171 


-183 


-189 




3.0 


0.657 X 10"^ 


-0.890 X 10-' 


0.950 X 10-' 


0.539 X 10 


-176 


-141 


-12 


-11 


-82 


-120 


-141 


-156 


-166 


-169 




3.5 


0.260 X 10-^ 


-0.138 X 10-' 


0.343 X 10-5 


A .5 1 A 1 A— 7 

—0.319 X 10 


-3014 


-2614 


-735 


-79 


-78 


-76 


-78 


-80 


-82 


-83 


4 


1.0 


0.836 X 10"^ 


0.251 X 10-2 


-0.670 X 10-" 


A 1 01 1 A— 2 


-39 


-40 


-48 


-74 


-112 


-144 


-159 


-166 


-175 


-180 




2.0 


0.820 X 10-^ 


0.205 X 10-2 


-0.872 X 10-' 


0.628 X 10-' 


-1 


-6 


-34 


-78 


-117 


-144 


-157 


-168 


-178 


-183 




3.0 


0.774 X 10"^ 


0.125 X 10-2 


-0.101 X 10-2 


0.169 X 10-' 


102 


84 


-1 


-79 


-121 


-144 


-158 


-171 


-182 


-186 




4.0 


0.708 X 10"^ 


0.197 X 10-' 


-0.545 X ur' 


0.780 X nr" 


330 


304 


105 


-67 


-122 


-136 


-150 


-163 


-174 


-177 




5.0 


0.518 X 10~^ 


-0.515 X IQ-' 


0.240 X 10-" 


-0.414 X 10-" 


-3170 


-2691 


-481 


-53 


-90 


-103 


-115 


-126 


-132 


-135 




5.5 


0.173 X 10-2 


-0.236 X 10-" 


0.145 X 10-* 


-0.323 X 10-' 


-9898 


-8248 


-1336 


-187 


-175 


-173 


-174 


-175 


-177 


-177 


6 


1.0 


0.836 X 10-2 


0.270 X 10-2 


-0.959 X 10-" 


0.142 X 10-2 


-41 


-41 


-48 


-74 


-113 


-145 


-159 


-166 


-174 


-179 




3.0 


0.805 X 10-2 


0.175 X 10-2 


-0.837 X 10-' 


0.176 X 10-' 


41 


29 


-27 


-82 


-119 


-143 


-158 


-169 


-178 


-182 




4.0 


0.765 X 10-2 


0.128 X 10-2 


-0.105 X 10-2 


0.134 X 10-' 


181 


151 


10 


-84 


-125 


-146 


-159 


-171 


-181 


-186 




5.0 


0.736 X 10-2 


0.417 X 10-' 


-0.650 X 10-' 


0.889 X 10 " 


433 


384 


97 


-79 


-126 


-138 


-152 


-165 


-175 


-179 




6.0 


0.676 X 10-2 


-0.873 X IQ-' 


0.531 X 10-" 


-0.115 X 10-5 


-1319 


-1044 


-112 


-63 


-97 


-115 


-132 


-147 


-157 


-160 




7.0 


0.258 X 10-2 


-0.515 X 10-" 


0.467 X 10-*' 


-0.153x10-* 


-10115 


-8462 


-852 


-213 


-209 


-212 


-214 


-219 


-222 


-223 


8 


1.0 


0.840 X 10-2 


0.276 X 10-2 


-0.536 X 10-" 


0.145 x 10-2 


-42 


-42 


-49 


-73 


-113 


-145 


-160 


-166 


-174 


-179 




3.0 


0.818 X 10-2 


0.200 X 10-2 


-0.749 X 10-' 


0.346 x 10-' 


14 


5 


-37 


-82 


-118 


-142 


-157 


-168 


-177 


-181 




5.0 


0.768 X 10-2 


0.120 X 10-2 


-0.104 X 10-2 


0.136 X 10-' 


267 


218 


20 


-87 


-126 


-146 


-159 


-171 


-181 


-186 




6.0 


0.757 X 10-2 


0.398 X 10-' 


-0.615 X 10-' 


0.797 X 10-" 


499 


441 


88 


-84 


-126 


-138 


-152 


-165 


-175 


-179 




7.0 


0.725 X 10-2 


-0.891 X 10-' 


0.435 X 10-" 


-0.301 X lO-" 


-893 


-679 


-40 


-73 


-101 


-120 


-138 


-153 


-164 


-167 




8.5 


0.342 X 10-2 


-0.704 X 10-" 


0.645 X lO-"" 


-0.211 X 10-* 


-5137 


-12194 


-990 


-271 


-273 


-278 


-284 


-290 


-294 


-295 
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